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Section  1 


PROGRAM  SWMAR7  FOR  TBE 
KSNIf  TAPERED-WIGGLER  PEL  OSCILLATOR  PROGRAM 


1.1  OVERVIEW 

ttie  basic  goal  of  this  program  is  to  provide  data  and  analysis  for 
assessing  the  potential  of  the  tapered-wiggler^ 1-1 )  FEL  oscillator  as  a 
high-efficiency  source  of  coherent  radiation.  The  advantage  of  the 
tapered -wiggler  FEL,  as  compared  to  the  original^ 1-2 )  fixed  pitch  device 
demonstrated  at  Stanford,  is  that  of  high  electron  kinetic-energy 
extraction  in  a  single  pass.  It  is  expected  that  the  tape red -wiggler 
concept  will  lead  to  development  of  high-efficiency  free-electron  lasers  if 
an  e-beam  recovery  stage  is  added  downstream  of  the  wiggler. 

The  fact  that  large  deceleration  of  electrons  could  be  achieved  with 
a  tape red -wiggler  design  was  first  demonstrated  in  the  Mathematical 
Sciences  Northwest/Boeing  Aerospace  Company  program,  and  since  then 
MSNW/BAC  have  measured^ 1-3 ^  extraction  as  large  as  4  percent  in  an 
amplifier  configuration,  and  recently  1  percent  has  been  realized^ 1-4 ^  in 
oscillators.  The  amplifier  gain  is  too  low  to  be  of  practical  interest, 
hence  realistic  systems  will  probably  be  configured  as  high-extraction 
oscillators.  There  is  therefore  incentive  to  determine  whether  the 
oscillator  can  be  made  to  work  at  extraction  comparable  to  that  achieved  in 
amplifier  experiments.  This  AFOSR  program  is  directed  toward  bridging  the 
gap  between  the  current  single  pass  measurements  and  high-extraction 
oscillators,  it  embodies  complementary  experimental  and  theoretical  tasks 
which  are  intended  to  assess  oscillator  potential  in  the  near  term  using 
single-pass  data. 

The  experimental  task  involves  measurement  of  the  FEL  interaction  as 
a  function  of  several  basic  parameters  for  the  purpose  of  verifying 
understanding  of  the  basic  interaction.  These  parameterization  experiments 
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were  completed  for  high  optical  flux  conditions  with  measurements  of 
extraction  as  a  function  of  electron  beam  energy  and  extraction  as  a 
function  of  optical  power.  The  data  obtained  confirms  elements  of  models 
used  to  predict  single-pass  oscillator  performance . 

The  theoretical  work  conducted  under  this  program  has  provided  some 
of  the  first  insight  into  the  transverse-mode  structure  in  tapered-wiggler 
PEL  oscillators.  The  analysis  developed  has  now  been  applied  to  the  near- 
concentric  cavity  geometries  of  interest  for  high  average  power  FELs,  and 
also  extended  to  three  dimensions.  This  3-0  model  has  been  used  to  analyze 
cavity  alignment  tolerances  including  the  effects  of  diffraction.  An 
interesting  result  of  the  analysis  is  the  finding  that  near-concentric  PEL 
cavities  cam  have  alignment  tolerances  which  are  much  less  stringent  than 
predicted  by  geometrical  (ray  optics)  analyses.  This  is  particularly 
important  in  view  of  the  extremely  small  tolerance  of  the  near-concentric 
cavity . 

An  additional  model  hats  been  developed  to  study  longitudinal  mode 
structure.  This  model  is  used  to  study  a  predicted,  but  not  yet  observed, 
instability  known  as  the  sideband  instability/ 1-5  ^  This  instability  is 
important  in  that  it  threatens  to  limit  the  electron  energy  extraction  in 
tapered-wiggler  oscillators.  The  sideband  instability  results  from  a 
resonance  between  the  rocking  frequency  of  electrons  in  the  ponderomotive 
potential  well,  and  the  beat  wave  produced  from  the  optical  wave  and  its 
sidebands.  The  sideband  modulates  the  optical  wave  in  such  a  manner  as  to 
increase  the  amplitude  of  the  rocking  in  the  potential  well,  leading 
eventually  to  detrapping.  Results  of  this  sideband  analysis  have  recently 
been  obtained  for  parameters  appropriate  for  high-power,  visible-wavelength 
systems.  Simulations  show  sideband  growth  with  associated  loss  of 
extraction  by  one-half  or  more.  The  instability  can  be  suppressed, 
however,  by  use  of  frequency-selective  optical  elements  in  the  laser 
resonator,  with  recovery  of  full  extraction. 


Ttoe  remainder  of  Section  1  contains  a  brief  summary  of  the 
experimental  work,  its  interpretation,  and  the  transverse-mode  modeling. 
Details  of  these  topics  have  been  published,  and  for  convenience  these 
publications  are  included  in  this  report  as  appendicies.  Appendix  A 
describes  the  MSNW/BAC  single-pass  amplifier  experiment  used  for  the 
parameter izat ions.  The  model  developed  to  fit  the  parameter izat ions  is 
described  in  Appendix  B.  Appendix  C  contains  the  transverse-mode  analysis, 
and  the  use  of  that  analysis  to  identify  diffractive  effects  on  cavity 
alignment  tolerances  is  given  in  Appendix  D.  Publications  and 
presentations  resulting  from  research  conducted  under  this  contract  are 
listed  in  Appendix  E. 

The  longitudinal-mode  analysis  used  for  the  sideband  instability  work 
has  not  yet  been  published,  and  is  described  in  detail  in  Section  2  of  this 
report.  In  addition  to  the  modeling,  practical  considerations  for 
providing  frequency  selectivity  in  PEL  optical  cavities  are  addressed  in 
Section  3. 

Section  2  deals  with  development  of  the  basic  numerical  model  of  the 
instability  and  calculation  of  suppression  afforded* by  various  cavity 
schemes.  The  calculations  differ  from  others  in  that  the  cause  of  long 
length  wigglers  and  long  electron  micropulse  length  is  treated,  these  being 
the  most  probable  parameters  for  near-term  visible-wavelength  high- 
extraction  systems.  The  suppression  calculations  are  made  for  both  edge 
and  bandpass  filter  functions  of  various  slopes,  and  show  the  degree  of 
filtering  necessary  to  achieve  a  given  suppression. 

Section  3  is  a  discussion  of  wavelength-dispersive  optical  elements. 
One  straightforward  approach  considered  is  the  use  of  wavelength-selective 
dielectric  coatings  on  the  end  mirrors.  Another  promising  technique 
involves  insertion  of  a  birefringent  filter  into  the  cavity.  Spatially- 
dispersive  schemes  are  also  considered.  Ibis  includes  gratings  and  a  novel 
approach  based  on  optical  mode  size  variation  with  wavelength. 


The  experimental  tasks  involve  parameterization  of  the  PEL 
interaction  in  terms  of  parameters  useful  to  oscillator  calculations. 

Single  pass  measurements  of  the  variation  of  the  PEL  interaction,  with 
respect  to  photon  intensity  and  photon  wavelength,  have  been  made  at  high 
flux  levels  using  the  MSNW/BAC  10  amplifier  hardware . ^ 1-3  ^  The 
wavelength  parameterization  is  necessary  because  of  an  expected  tendency  of 
tapered-wiggler  devices  to  chirp  during  start-up,  and  it  has  been  measured 
on  a  single-pass  basis  not  by  variation  of  wavelength,  but  by  variation  of 
the  electron  energy.  The  interaction  can  be  measured  in  terms  of  either 
the  net  electron  energy  loss  or  equivalently  the  photon  gain.  Measurement 
of  the  energy  loss  has  been  made  over  a  range  of  photon  intensities 
spanning  the  onset  of  trapping  to  saturation. 

The  parameter izat ions  obtained  are  shown  in  Figures  1-1  and  1-2.  The 
electron  energy  extraction,  as  a  function  of  e-beam  energy,  is  shown  in 
Figure  1-1 .  Peak  extraction  is  observed  near  19  MeV.  When  the  electron 
energy  is  detuned  below  18.4  MeV,  net  acceleration  of  the  electrons  is 
seen.  Electron  deceleration  and  energy  extraction  is  observed  over  a  3 
percent  range  in  energy,  corresponding  to  a  6  percent  equivalent  range  in 
the  optical  spectrum.  The  data  are  in  good  agreement  with  the  theoretical 
curves  drawn,  under  the  assumption  of  perfect  focusing,  perfect  alignment, 
and  a  diffraction-limited  optical  beam.  Emittance  (defined  as  yirzB ) 
values  of  o  and  0.023  cm-rad  are  used.  The  latter  value  is  consistent  with 
the  estimated  combined  effect  of  emittance  and  misalignment.  The  energy 
scale  of  the  theoretical  curves  has  been  downshifted  1.5  percent  to  provide 
an  improved  fit.  This  shift  is  within  the  uncertainty  of  the  electron 
spectrograph  calibration.  The  points  shown  do  not  include  all  the  data 
taken.  With  misalignment,  poor  focusing,  or  other  problems,  it  is  always 
possible  to  achieve  results  in  Which  the  extraction  values  are  low,  but  it 
is  not  possible  to  achieve  extraction  results  that  are  artifically  high 
except  by  actual  measurement  errors. 
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1.  Extraction  as  a  Function  of  Input  Energy. 
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Figure  1-2 .  Extraction  as  a  Function  of  Laser  Power  in  the 
Wiggler. 
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Electron  energy  extraction  has  also  been  measured  as  a  function  of 
laser  power  as  shown  in  Figure  1-2.  The  extraction  is  seen  to  increase 
with  laser  power  as  expected.  The  onset  of  trapping  is  predicted  at 
roughly  50  Mr,  and  the  trapping  fraction  increases  steadily  to  about  60 
percent  at  500  MW.  The  data  shown  is  taken  from  three  time-resolved 
electron  energy  spectrum  records.  Since  the  laser  power  is  also  measured 
as  a  function  of  time,  it  is  a  simple  matter  to  determine  extraction  as  a 
function  of  power.  Again,  perfect  focusing  and  alignment  are  assumed  for 
the  theoretical  curve,  and  a  0. 023  cm-rad  normalised  emittance  is  used  as  a 
rough  estimate  of  combined  emittance  and  misalignment  effects.  When  the 
calculation  is  repeated  with  zero  emittance,  the  result  is  a  curve  parallel 
to  that  shown  and  about  1  percent  higher  extraction. 

1.3  TRANVERSE-MOOE  STRUCTURE  SUMART 


The  PEL  requires  good  transverse  optical  mode  quality  to  provide  a 
uniform,  high  photon  intensity  within  the  wiggler  and  thereby  maximize  the 
interaction  with  the  electrons.  In  addition,  production  of  a  nearly 
diffraction-limited  beam  will  imply  good  output  beam  quality.  Recent 
calculations  of  mode  evolution  in  injected  linear  cavities  are  sumnarized 
as  follows.  The  code  numerically  solves  Maxwell’s  equations  for 
cylindrically-symmetric  geometries.  It  allows  one  to  follow  development  of 
the  mode  structure  of  the  tape red -wiggler  oscillator,  starting  with  an 
initial  injected  wave  at  saturation  intensity.  This  injected  wave  develops 
over  many  round  trips  of  the  optical  cavity  according  to  the  influence  of 
trie  PEL  interaction,  diffraction,  and  interaction  with  other  elements  of 
the  optical  cavity.  The  e-beam  radius  is  typically  less  than  that  of  the 
photon  beam  and  the  size  mismatch  drives  higher  order  modes  in  addition  to 
te*Vk>-  Of  additional  interest  is  the  effect  of  apertures  on  the  mode 
structure,  especially  those  associated  with  the  wiggler  magnets.  E-field 
truncation  at  the  ends  of  the  wiggler  causes  mode-dependent  cavity  losses 
which  of  course  are  higher  for  higher-order  modes.  The  truncation  also 
causes  mode  mixing,  Which  transfers  power  from  low  to  high  order  modes. 
Optimization  studies  have  shown  that  the  minimum  wiggler  bore,  consistent 
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with  acceptable  distortion  due  to  clipping,  yields  the  largest  PEL 
interaction . 

Evolving  phase  and  amplitude  profiles  of  the  optical  beam  are 
calculated  by  direct  integration  of  Maxwell's  equations  with  the  electrons 
providing  the  driving  term.  A  convenient  tool  for  understanding  this 
evolution  is  the  projection  of  the  optical  wave  into  normal  modes  of  the 
optical  cavity.  That  is,  the  fraction  of  total  optical  power  in  any 
particular  cavity  mode  can  be  observed  as  a  function  of  the  round-trip 
number.  One  sees,  for  example,  that  a  pure  TEMqq  wave  injected  into  the 
optical  cavity  at  the  start  of  the  calculation  evolves  an  appreciable 
fraction  of  TEM10  mode  over  the  course  of  several  round  trips.  This  mode 
mixing  is  caused  primarily  by  the  nonuniform  gain  medium,  and  to  a  lesser 
extent,  by  the  presence  of  the  wiggler  entrance  and  exit  apertures. 

The  mode  evolution  during  the  first  35  round  trips  after  TEMqq  mode 
injection  is  shown  in  Figure  1-3  for  three  different  cavity  lengths,  each 
with  10  percent  output  coupling.  Higher-order  mode  content  is  especially 
evident  in  confocal  cavities  (those  with  mirror  separation  equal  to  the 
radius  of  curvature),  and  concentric  cavities  (those  with  mirror  separation 
equal  to  twice  the  radius  of  curvature ) .  These  cavities  can  support 
unusual  mode  structure  because  the  relative  phase  slippage  between  cavity 
modes  over  one  round  trip  is  an  integral  multiple  of  2V ,  allowing 
constructive  interference  between  higher-order  modes  produced  on  each  round 
trip.  As  the  cavity  length  is  changed  away  from  the  confocal  or  concentric 
condition  without  changing  the  degree  of  mode  selectivity  due  to 
aperturing,  the  fraction  of  mode  decreases  dramatically. 

Significant  higher  order  mode  content  may  also  be  found  in  high-gain 
systems.  As  shown  in  Figure  1-4,  higher  gain  systems  will  tend  to  exhibit 
unusual  mode  structure  even  at  intermediate  cavity  lengths.  Systems  with 
high  steady-state  gain  also  have  large  output  coupling  and  therefore  have  a 
short  ring  down  time.  Constructive  interference  between  higher-order  modes 
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produced  over  many  round  tripe  ie  not  required  for  significant  higher-order 
■ode  content  in  these  systems. 

Figure  1-5  cospares  the  steady-state  on-axis  intracavity  photon 
intensity  with  that  of  the  initially  injected  TEHqq  mode  in  a  near- 
concentric  cavity.  Bigh  average  power  FELs  with  linear  cavities  require 
near-concentric  cavities  due  to  mirror  loading.  The  rapid  variations  in 
intensity  are  due  to  diffraction  effects  from  the  truncation  at  the  ends  of 
the  wiggler.  The  truncation  occurs  at  radius  h  *-  1.8  w,  where  w  is  the  1/e 
point  in  E-field  for  a  TEMq 0  wave  at  the  aperture.  The  TEM10  mode  content 
results  in  the  striking  asymmetry  between  the  forward  and  backward  moving 
waves  in  the  wiggler.  Such  an  asymmetry  can  be  supported  in  the  FEL  since 
the  gain  mechanism  is  active  in  one  direction  only. 

The  quality  of  the  output  beam  is  excellent  for  each  of  the  cases 
studied,  being  of  nearly  diffraction-limited  quality.  This  result  is 
somewhat  surprising  in  view  of  the  TEM10  mode  content  which  complicates  the 
intracavity  structure,  but  in  reality  the  higher-mode  content  is  an 
indication  of  mode  mismatch  within  the  cavity,  rather  than  wave  front 
aberration.  For  the  example  shown  in  Figure  5,  the  focused  output  beam  has 
a  Strehl  ratio  of  about  98  percent. 

A  three-dimensional  version  of  this  wave  front  propagation  analysis 
has  been  developed  and  can  be  used  for  mode  analysis  similar  to  that  shown 
previously,  but  additional  important  effects  such  as  misalignment  earn  be 
included.  In  addition,  the  3-0  code  provides  for  analysis  of  complex 
cavities  such  as  ring  cavities  employing  glancing- incidence  mirrors 
intended  for  use  at  high  average  power.  Results  of  including  diffraction 
in  mirror  alignment  tolerance  analysis  of  linear  cavities  are  given  in 
Appendix  0.  Tolerance  to  cavity  misalignment  is  generally  studied  with 
geometrical  optics  codes,  but  the  low  Fresnel  number  FEL  cavities  are 
dominated  by  diffraction  and  one  finds  that  the  alignment  tolerances  do  not 
follow  conventional  rules  based  on  ray  tracing. 
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Section  2 


SIDEBAND  SUPPRESSION  SIMULATION 

A  simulation  model  has  been  developed  for  study  of  the  potentially 
serious  Raman  sideband  instability  first  described  by  Kroll  and 
Rosenbluth.^ 2-1 ^  The  instability  is  characterised  by  development  of 
longitudinal  amplitude  and  phase  modulation  in  the  optical  pulse  due  to  the 
generation  of  new  frequency  components  or  ’’sidebands"  in  the  laser 
spectrum.  The  instability  is  expected  to  become  evident  in  FEL  oscillators 
at  optical  powers  high  enough  to  trap  electrons  in  the  ponderomotive 
potential  well  of  the  FEL  interaction,  because  synchrotron  oscillations  of 
trapped  electrons  lead  to  axial  modulation  of  the  gain. 

The  instability  is  predicted  to  result  in  a  loss  of  electron  trapping 
efficiency  in  multi-pass  FELs  using  highly- tape red  wigglers.  Simulations 
of  untapered  wigglers^ 2-2 ' 2-3 ^  show  that  the  sideband  instability  actually 
leads  to  enhanced  extraction  by  promoting  chirp  to  the  frequency  of  peak 
saturated  gain.  But  tapered  wigglers  rely  upon  trapping  of  electrons  in 
decelerating  ponderomotive  potential  wells  to  achieve  enhanced  extraction 
efficiency.  Highly-tapered  wigglers  trap  the  electrons  for  many 
synchrotron  periods.  The  sideband  instability  can  cause  detrapping  within 
a  distance  as  short  as  one  synchrotron  period,  thus  potentially  leading  to 
severe  loss  of  extraction  in  highly-tapered  systems. 

This  study  has  provided  the  first  simulation  of  the  time-dependent 
sideband  evolution  for  parameters  of  a  visible  long-pulse  tapered-wiggler 
FEL  oscillator,  where  there  are  three  or  more  synchrotron  periods  within 
the  wiggler  length.  Simulations  show  sideband  growth  with  associated  loss 
of  extraction  efficiency  by  one-half  or  more.  The  same  simulation  shows 
that  the  instability  can  be  suppressed  by  use  of  frequency-selective 
elements  in  the  optical  cavity,  with  full  recovery  of  the  extraction 
efficiency.  These  analyses  will  be  highly  relevant  for  interpretation  of 


experimental  results  as  they  become  available  and  for  guidance  of  future 
oscillator  design  work. 


2.1  DESCRIPTION  OP  MODEL 


the  mathematical  basis  for  the  sideband  instability  has  been 
summarized  by  Goldstein  and  Colson. ^ 2-2 ^  the  electron  motion  in  the 
ponderomot i ve  well,  for  small  deviations  from  the  resonant  phase  and  an 
unmodulated  E-field,  may  be  described  by  a  harmonic  oscillator  equation, 
this  implies  that  such  electrons  might  be  expected  to  couple  to  light  whose 
wavelength  \'B  is  slightly  shifted  from  the  original  resonant  wavelength 


[2-1] 


where  Xw  is  the  wiggler  wavelength.  is  the  spatial  period  of  the 
harmonic  behavior  for  an  electron  near  the  bottom  of  the  bucket 
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G  «  J0(y)-Jj(y)  is  the  coupling  factor  for  a  planar  wiggler  given  by  the 

2  2 

difference  between  two  Bessel  function  factors,  y  *  a^/2(l+aw),  ^ r  is  the 
synchronous  phase  angle  for  trapped  electrons 


and  the  optical  electric  field  amplitude  eq  and  peak  wiggler  magnetic  field 
B0  are  measured  by  (cgs  units) 
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respectively.  The  quantity  e  is  the  minimum  E-field  for  electron  trapping 
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for  a  wiggler  of  resonant  energy  yT,  fractional  energy  taper  A>r/yr,  and 
length  1^.  In  actuality,  many  electrons  undergo  large  amplitude  phase 
oscillations.  Which  means  that  their  motion  is  described  by  a  more  complex 
nonlinear  pendulum  equation.  Such  electrons  couple  to  a  continuum  of  light 
waves  of  different  frequencies  close  to  that  given  by  Equation  [2-1]. 
Kroll's  analysis^ 2-1  *  predicts  a  spectrum  of  unstable  waves  with  the 
largest  growth  rates  for  waves  satisfying  Equation  [2-1]  with  the  positive 
sign.  Thus  one  expects  gain  for  lower  frequency  sidebands  (longer 
wavelengths)  and  absorption  for  higher  frequencies.  Another  way  to 
describe  this  phenomenon  is  that  light  at  X8  is  Raman  shifted  to  Xg  >  Xg 
along  with  excitation  of  increased  sloshing  of  electrons  in  the 
ponderomotive  well. 

The  spatial  modulation  period,  Xm,  associated  with  generation  of  a 
sideband  satisfying  Equation  [2-1 J  is  given  by 
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This  is  to  be  compared  with  the  slippage  distance 
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Which  is  the  distance  by  which  a  resonant  electron  slips  back  relative  to  a 
plane  wave  in  one  transit  through  the  wiggler.  Most  previous  simulations 
of  sideband  evolution^ 2-4 ^  have  observed  initial  onset  of  modulation  with  a 
period  approximately  equal  to  the  slippage  distance,  which  corresponds  to  a 
synchrotron  period  comparable  to  the  wiggler  length.  Evidently  sideband 
growth  does  not  occur  until  the  optical  power  grows  to  the  point  where  the 
synchrotron  period  becomes  comparable  to  or  shorter  than  the  wiggler 
length.  The  fractional  frequency  shift  of  the  sideband  is  then  of  order 
1/N,  where  W  is  the  number  of  wiggler  periods.  These  simulation  studies 
have  considered  short,  mi Idly- tapered  wigglers  in  which  there  is  only  about 
one  synchrotron  period  within  the  wig-ler  length. 


Linearly-tapered  wigglers  may  be  characterized  by  the  single 
dimensionless  parameter^ 2-5 ' 2-6 ^ 
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This  work  differs  from  others  primarily  in  that  we  consider  a  long,  highly- 
tapered  wiggler  appropriate  for  a  high-efficiency  visible-wavelength  FEL, 
for  which  6  »  12577 .  This  is  an  order  of  magnitude  larger  than  the  mildly- 
j  tapered  wigglers  of  0  <  1577  for  which  simulations  of  time-dependent 

oscillation  behavior  for  long  electron  pulses  have  previously  been 
presented  in  the  literature. ^ 2-7 ' 2-8 ^  Examination  of  Equations  [2-2],  [2- 
3],  and  [2-6]  reveal3  that  the  number  of  synchrotron  periods  within  the 
wiggler  length  is  given  by 
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Thus  we  see  that  at  any  given  resonant  phase  angle  (fT  *  40  degrees 
maximizes  the  product  of  electron  trapping  fraction  and  deceleration  rate 
per  unit  E-field^ 2-9 ^ ),  the  wiggler  considered  in  this  study  has  roughly 
three  times  as  many  synchrotron  periods  within  the  wiggler  length.  For 
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such  long,  highly-tapered  wigglers,  the  fractional  frequency  shift  of  the 
sideband  given  in  Equation  [2-1]  will  be  approximately 


[2-11] 


Time-dependent  sideband  evolution  in  longer,  more  highly- tapered 
wigglers  has  previously  been  simulated/2-10)  but  only  under  the  assumption 
of  an  electron  pulse  length  comparable  to  or  shorter  than  the  slip  length. 
The  codes  used  may  be  applied  to  the  much  longer  pulse  experiments  of 
interest  to  DOD  goals  in  which  the  pulse  length  is  20  or  more  slippage 
distances/2-7^  but  the  computer  time  requirements  are  very  large.  In  this 
work,  we  consider  very  long  pulses  by  use  of  periodic  boundary 
conditions^2-3'2-8 )  to  examine  the  pulse  modulation.  The  periodic  boundary 
condition  model  allows  treatment  of  a  short  section  of  much  longer  electron 
and  optical  pulses,  thus  greatly  reducing  computer  time  requirements. 

The  modulation  of  the  optical  pulse  is  studied  directly  in  the  space- 
time  domain.  The  one-dimensional  model  for  the  time-dependent  evolution  of 
the  optical  field  modulation  is  developed  under  the  assumptions  of  a  plane- 
polarized  tapered  wiggler  field 
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and  diffaction-limited  plane-polarized  optical  field 
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Arbitrary  wiggler  tapering  is  provided  by  prescribing  the  axial  variation 
in  the  magnetic  field  amplitude  BQ.  The  spatial  variation  in  amplitude  and 
phase  for  a  TEMq0  beam  of  Rayleigh  range  ZR  focused  at  the  wiggler  center 
is  given  by  the  functions^ 2_ 11 ) 
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where  q  »  1^/27^,  z  is  measured  from  the  wiggler  entrance,  the  1/e 
amplitude  radius  of  the  beam  is 
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Optical  field  evolution  due  to  the  PEL  interaction  is  included  by  the  axial 
variation  in  the  E- fie Id  amplitude  eso  and  phase  0  explicitly  shown  in  the 
equations. 

The  electron  beam  is  assumed  to  also  be  axi symmetric  and  to  have  a 
parabolic  density  profile 


[2-18] 


where  I  is  the  beam  current  and  reb  is  the  beam  radius.  The  parabolic 
density  profile  is  a  close  approximation  to  the  distribution  resulting  from 
uniformly-filled  emittance  phase  space.  The  beam  radius  is  related  to  the 
normalized  emittance  eN  --  yvxQ  by 
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is  the  wavenuihber  for  betatron  oscillations  resulting  from  distributed  two- 
plane  focusing^ 2-12  )  in  the  wiggler. 


Electron  energy  spread  and  emittance  axe  important  effects  Which  tend 
to  reduce  the  PEL  interaction  strength  for  the  visible  PEL  parameters 
considered  here. ^ 2-12 )  Energy  spread  cam  be  included  directly  in  the  model 
by  introducing  a  dispersion  of  the  initial  energies  of  the  sample  electrons 
used  to  drive  the  optical  wave.  However,  it  is  not  possible  to  directly 
include  transverse  beam  effects  such  as  emittance  in  a  1-0  axial  model. 
E-beam  emittance  results  in  loss  of  interaction  strength  due  to 
introduction  of  both  effective  electron  energy  spread  and  less  than  ideal 
overlap  of  the  electron  and  photon  beams.  For  the  conditions  considered 
here,  loss  of  overlap  is  of  far  more  significance  than  effective  energy 
spread . ^ 2-12 )  Consequently,  effective  energy  spread  has  been  neglected  but 
overlap  effects  have  been  included  in  an  approximate  way  using  appropriate 
weighted  averages.  The  column-averaged  E-field  amplitude  and  phase  seen  by 
the  e-beam  are  given  by 


» 


9 


9 


I 


I 


2-7 


I 


f  *Q(r,z)  ne(r,z) 


2ffrdr 


[2-22] 


J  ^  r )  2»rdr 
J  o 


» (-4H  -“"'ItH 


In  the  limit  of  zero  e-beam  radius,  ega  reduces  to  eB0w0/w .  The  Lorentz 
force  equation  describing  the  motion  of  a  sample  electron  in  the 
ponderomotive  potential  well  formed  by  the  wiggler  and  optical  fields  may 
be  written  in  terms  of  the  column-averaged  quantities 
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where  t  «■  (kw  kg)z  -  wgt  is  the  phase  of  the  electron  relative  to  a 
plane  wave. 

The  evolution  of  the  optical  wave  is  determined  based  on  the  self- 
consistent  interaction  with  a  number  of  sample  electrons.  The  sample 
electrons  are  injected  into  the  wiggler  on  each  pass  uniformly  distributed 
in  f  and  distributed  in  y  according  to  the  energy  spread  of  the  electron 
beam.  In  the  slowly  varying  phase  and  amplitude  approximation/2-13^  the 
driven  wave  equation,  in  a  reference  frame  moving  with  the  photons,  reduces 
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where  f  is  a  gsostetrical  "fill  factor" 
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and  the  averages  indicated  by  the  angular  brackets  are  made  over  the  set  of 
sample  electrons. 

We  study  the  evolution  of  a  long  optical  pulse  by  considering  a 
length  W  near  the  center  of  the  pu]se  and  imposing  periodic  boundary 
conditions  at  each  end.  The  window  width  is  chosen  to  be  at  least  as  long 
as  the  slip  length,  since  this  is  an  important  characteristic  distance  in 
the  problem.  Because  of  the  periodic  boundary  conditions,  the  optical 
pulse  may  be  decomposed  into  a  set  of  discrete  optical  frequencies.  The 
choice  of  the  width  of  the  window  between  the  periodic  boundaries  and  the 
number  of  spatial  grid  points  within  the  window,  Np,  defines  the  discrete 
frequencies  which  are  handled  in  the  model. 
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The  center  wavelength  is  taken  to  be  the  resonant  wavelength,  defined  to  be 
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where  y  is  the  mean  initial  electron  energy  and  a^,  is  the  aw  value  at 
the  wiggler  entrance.  For  a  window  width  of  one  slip  length  and  a  system 
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with  four  synchrotron  periods  along  the  wiggler  length,  the  sideband  may  be 
expected  to  be  separated  from  the  main  line  by  four  of  the  discrete 
frequency  intervals.  Longer  windows  provide  better  frequency  resolution  at 
the  expense  of  increased  computer  time  requirements. 

The  electron  slippage  is  included  by  allowing  the  sample  electrons  to 
drift  across  the  optical  field  window  as  the  electron  and  optical  pulses 
propagate  down  the  wiggler.  In  accordance  with  the  periodic  boundary 
condition  model,  those  electrons  which  slip  out  of  one  side  of  the  window 
reappear  at  the  opposite  side  of  the  window.  In  order  to  simplify  the 
bookkeeping,  the  rate  of  slippage  for  all  electrons  is  taken  to  be  the  rate 
for  a  resonant  electron.  This  approximation  neglects  the  actual 
distribution  of  electron  axial  velocities,  which  may  be  fairly  appreciable 
in  a  tapered  wiggler.  The  difference  in  slippage  lengths  between  trapped 
and  un trapped  electrons,  for  example,  is 


For  the  cases  shown  here,  the  taper  is  about  10  percent.  Nevertheless, 
this  approximation  appears  to  be  acceptable  considering  that  untrapped 
particles  appear  to  play  a  very  minor  role  in  the  sideband  instability. 

Evolution  of  the  modulation  due  to  the  sideband  instability  is 
studied  by  solving  the  driven  wave  equation  within  the  window  over  many 
passes  through  the  oscillator.  After  each  pass  the  E- field  amplitude  is 
renormalized  to  account  for  round-trip  losses  such  as  output  coupling. 
After  many  passes  the  optical  pulse  reaches  a  steady  state  in  which  the 
8 ingle- pas s  gain  equals  the  losses. 

2.2  SIMULATION  OF  SIDEBAND  EVOLUTION 


Oscillators  incorporating  long,  highly-tapered  wigglers  can  be 
expected  to  be  subject  to  somewhat  more  severe  detrapping  should  sidebands 


appear,  since  there  is  more  opportunity  for  the  instability  to  drive 
sloshing  in  the  potential  well.  Figure  2-1  presents  results  of  an  early 
study  to  evaluate  the  possible  extraction  degradation,  should  sidebands 
appear.  This  calculation  is  not  self-consistent  in  that  it  does  not 
predict  the  magnitude  or  wavelength  offset  of  the  sidebands.  It  is  simply 
assumed  that  a  single,  well-separated  sideband  is  present  at  a  very  modest 
power  level,  4  percent  of  the  power  at  the  carrier  frequency.  The  E-field 
of  the  sideband  is  then  20  percent  of  the  carrier  E-field.  Figure  2-1 
shows  the  electron  energy  extraction  as  a  function  of  the  wavelength  offset 
of  the  sideband,  calculated  for  a  single  pass  through  each  of  two  wigglers 
of  different  length.  For  the  longer  wiggler,  a  very  distinct  resonance  is 
found,  which  occurs  for  a  sideband  wavelength  offset  of  approximately 


This  offset  is  at  a  frequency  which  causes  the  modulation  due  to  the 
sideband  to  be  in  phase  with  the  synchrotron  oscillations  of  trapped 
electrons  and  is  in  agreement  with  the  predicted  sideband  offset  given  in 
Equation  [2-9].  The  instability  has  a  finite  bandwidth  because  of  the 
range  of  synchrotron  periods  due  to  anharmonicity  of  the  bucket  and  the 
range  of  aw  and  es  values  within  the  wiggler.  E-field  phase  modulation  due 
to  the  sideband  results  in  side-to-side  motion  of  the  ponderomotive  well. 

In  the  longer  wiggler,  the  shaking  of  the  well  induces  large  amplitude 
sloshing  of  trapped  electrons,  causing  detrapping  and  the  loss  of 
extraction  shown .  The  parameters  assumed  for  the  longer  wiggler  in  this 
calculation  are  very  similar  to  those  of  the  preliminary  design  of  a  high- 
extraction  visible  oscillator  experiment  described  in  Table  2-1. 

These  calculations  give  am  indication  of  the  possible  serious  impact 
of  the  sideband  instability  for  conditions  of  a  visible-wavelength  tapered- 
wiggler  PEL .  With  as  little  as  4  percent  of  the  total  optical  power  in  the 
sideband,  the  extraction  could  be  lowered  from  the  design  value  of  5 
percent  to  a  value  of  only  3  percent.  But  these  calculations  are  not  self- 
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Table  2-1 

PRELIMINARY  POINT  DESIGN  FOR  VISIBZ2 
OSCILLATOR  EXPERIMENT 


Electron  Beam 

Initial  Energy 
Peak  Current 
Pulse  Length 


120  MeV 
>100  A 
25  ps 


wriggler 

Length 

Taper 

Wavelength 

Peak  Field  Strength 

a 

w 

Optical  Cavity 

Resonant  Wavelength 
Rayleigh  Range 
Length 
Outcoupling 

Derived  Quantities*®* 

Slippage  Lengths  per  Pulse 

E-field  for  5  percent  Extraction,  eg 

Photon  Power  for  5  Percent  Extraction 

Synchrotron  Period  at  5  Percent  Extraction 

Expected  Number  of  Modulation  Periods 
per  Slip  Distance 

Expected  Frequency  Offset  of  Sideband 


5  m 

12  Percent 
2.02  cm 
10.0  T 
1.33 


0.5/un 
2.4  m 
60  m 
Variable 


65 

20  cm"  1 
3.3  GW 
1.34  m 
3.7 


1.5  Percent 


(a), 


Without  allowance  for  emittance  and  energy  spread. 


2-13 


consistent  in  that  they  do  not  include  the  influence  of  the  electron 
dynamics  on  the  optical  E-field  evolution. 

Results  of  a  calculation  of  the  self-consistent  evolution  of  the 
optical  spectrum  for  the  parameters  of  Table  2-1  are  shown  in  Figure  2-2. 

In  this  simulation  the  electron  beam  is  assumed  to  be  of  ideal  quality, 
that  is,  energy  spread  and  emittance  acre  neglected.  The  simulation  assumes 
an  electron  beam  current  of  200  A  and  round-trip  cavity  losses  (output 
coupling)  of  10  percent.  Wie  wiggler  taper  consists  of  a  0.75  m  uniform 
section  followed  by  a  4.25  m  linearly-tapered  section  of  12.5  percent 
resonant  energy  change.  The  inclusion  of  a  short  constant  section  in  the 
wiggler  prescription  allows  enhancement  of  the  small-signal  gain  to  above 
the  saturated  gain  level  at  the  expense  of  a  modest  reduction  in  small- 
signal  linewidth.  The  power  levels  shown  on  the  figure  refer  to  the 
instantaneous  power  of  the  optical  pulse  within  the  cavity.  The  initial 
optical  field  is  arbitrary.  To  roughly  approximate  the  initial  incoherent 
properties  of  the  E- field,  the  simulation  is  initiated  with  all  possible 
frequencies  present,  but  randomly  phased  with  respect  to  each  other.  The 
simulation  is  seeded  at  a  power  level  of  approximately  1  watt  per  frequency 
channel,  which  is  representative  of  the  spontaneous  emission  power  level. 
Use  of  various  representations  of  the  initial  spectrum  result  in  somewhat 
different  details  in  the  spectral  evolution  but  do  not  change  the 
qualitative  results.  A  window  width  of  twice  the  slippage  distance  and  64 
spatial  grid  points  are  used  in  the  simulation.  The  simulation  is  thus 
conducted  within  a  frequency  bandwidth  of  approximately  12  percent  full 
width . 

After  a  number  of  passes  through  the  oscillator,  the  laser  picXs  a 
narrow  line  from  the  initial  seed.  The  frequency  and  shape  of  the  line  is 
consistent  with  the  small-signal  gain  curve,  as  shown  in  Figure  2-2.  The 
power  level  of  this  line  exponentiates  at  the  small-signal  gain  rate  of 
about  65  percent  per  pass.  As  the  power  level  approaches  saturation,  the 
line  Chirps  slightly  by  growth  of  the  longer  wavelength  wing  of  the  line. 
Upon  reaching  a  power  level  of  several  GW,  sufficient  to  trap  electrons,  an 
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upper  sideband  begins  to  form.  By  pass  number  120,  shown  in  Figure  2-2, 
the  sideband  has  grown  to  the  point  where  the  electron  trapping  efficiency 
is  seriously  degraded.  The  sideband  offset  of  about  1.5  percent  is 
consistent  with  the  number  of  synchrotron  periods  ("4)  at  this  power 
level.  As  predicted  by  Kroll, ^ 2-1 ^  the  upper  sideband  is  dominant, 
although  both  upper  and  lower  spectral  features  appear.  The  E-field 
amplitude  and  phase  modulation  due  to  the  sideband  is  shown  in  Figure  2-3. 
The  modulation  period  is  about  one-fourth  the  slippage  distance  which  is 
consistent  with  Equation  [2-7].  The  sidebands  continue  to  evolve  to  higher 
wavelength  and  the  spectrum  develops  multiple  features. 

The  buildup  of  the  sideband  results  in  a  loss  of  trapping  efficiency 
and  a  corresponding  decay  of  the  laser  power.  As  shown  in  Figure  2-4,  a 
true  steady  state  is  not  reached,  lhe  laser  reaches  an  oscillatory  quasi¬ 
steady  state  in  which  it  vacillates  around  an  average  power  of  about  7  GW. 
The  power  oscillations  correspond  to  slow  growth  and  decay  of  various 
spectral  features.  The  extraction  efficiency  as  a  function  of  time  is 
shown  in  Figure  2-5.  The  loss  of  extraction  upon  buildup  of  the  sideband 
instability  is  clearly  shown  in  the  figure.  The  extraction  efficiency  in 
quasi-steady  state  averages  about  3  percent,  considerably  lower  than  the 
ideal  value  of  about  7  percent  which  would  be  obtained  if  the  laser  were 
operating  with  a  narrow  line  at  the  frequency  of  peak  gain.  As  may  be  seen 
by  comparison  with  the  saturated  gain  curve  for  a  single  narrow  line  of 
equal  power  shown  in  Figure  2-2,  the  operating  spectrum  is  not  only  complex 
but  at  a  longer  than  optimum  wavelength. 

Results  of  a  study  of  the  degradation  of  extraction  efficiency  due  to 
the  sidebands  as  a  function  of  output  coupling  is  shown  in  Figure  2-6. 

There  are  a  range  of  output  couplings  for  which  the  laser  will  evolve  to 
power  levels  sufficient  for  trapping.  The  ideal  extraction  efficiency  at 
saturation  for  operation  with  a  single  narrow  line  which  is  allowed  to 
chirp  is  indicated  by  the  solid  line.  When  the  effect  of  the  sidebands  is 
included,  the  extraction  iB  reduced  to  approximately  40  percent  of  the 
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Figure  2-3.  Optical  Field  Amplitude  and  Phase  at  Pass  Number  120 
Showing  Modulation  Due  to  Sideband  Instability. 
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expected  value.  Of  course,  sidebands  do  not  form  for  output  couplings 
which  are  too  high  to  allow  the  laser  to  achieve  trapping. 

2.3  S  INDUCTION  OF  SIDEBAND  SUPPRESSION  TECHNIQUES 

The  model  of  sideband  evolution  allows  study  of  sideband  suppression 
techniques.  Kroll  originally  suggested^ 2-1 ^  that  the  addition  of  electron 
energy  spread  would  reduce  the  sideband  growth  rate  since  the  growth  rate 
depends  on  the  details  of  the  phase  space  distribution  of  electrons  trapped 
in  the  potential  well.  Simulations  have  indeed  shown ^ 2-2 ^  that  large 
initial  energy  spread  reduces  the  growth  rate  but  does  not  eliminate  the 
instability.  Here  we  investigate  the  possibility  of  totally  eliminating 
sideband  activity  by  use  of  a  dispersive  element  in  the  optical  cavity  to 
provide  wavelength  selectivity. 

A  number  of  possible  schemes  for  provision  of  wavelength  selectivity 
are  considered  in  Section  3.  As  an  example  of  the  performance  of  these 
techniques,  the  simulation  shown  in  Figures  2-2  through  2-5  has  been 
repeated  including  optical  filters  in  the  cavity.  Figure  2-7  shows  the 
filter  functions  assumed  in  the  two  example  calculations.  The  round-trip 
cavity  loss  on  the  main  line  is  retained  at  10  percent.  In  Case  (a)  of 
Figure  2-7,  a  4  percent  bandpass  filter  is  placed  in  the  cavity.  We  define 
the  bandwidth  to  be  the  width  at  the  50  percent  loss  points.  The  sin3(Xs) 
dependence  of  the  cavity  loss  is  characteristic  of  that  produced  by  a 
birefringent  filter.  In  Case  (b),  an  edge  filter  iB  used  with  a  cutoff  at 
longer  wavelengths.  The  edge  filter  function  corresponds  to  the  wavelength 
selectivity  possible  by  operation  of  the  laser  near  the  edge  of  a  10 
percent  bandpass  filter  produced  by  a  multilayer  dielectric  mirror 
coating.  In  both  cases  the  selectivity  functions  are  configured  in  such  a 
way  that  the  losses  are  large  at  the  sideband  wavelength  offset  of  1.5 
percent  (4  to  5  times  the  loss  on  the  main  line). 

Results  of  these  simulations  are  compared  in  Figure  2-8.  Both 
filters  successfully  eliminate  the  sideband  instability.  The  optical 
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spectra  at  pass  number  300  are  very  narrow  compared  to  the  case  with  no 
frequency  selectivity.  Furthermore,  the  operating  wavelength  is 
approximately  the  wavelength  of  peak  saturated  gain,  rather  than 
overshooting  to  longer  wavelengths.  Consequently,  the  laser  saturates  at  a 
muCh  higher  power  level  and  the  extraction  efficiency,  as  shown  in  Figure 
2-9,  evolves  to  over  7  percent,  representing  very  efficient  electron 
trapping.  In  the  case  of  the  edge  filter,  the  temporary  drops  in 
extraction  shewn  in  the  figure  correspond  to  chirp  of  the  laser  spectrum 
from  one  discrete  frequency  to  another.  This  temporary  extraction  loss  is 
believed  to  be  a  nonphysical  result  due  to  the  limited  frequency  resolution 
of  the  model. 

An  additional  sisiulation  was  also  conducted  using  a  birefringent 
filter  function  with  a  full  width  of  6  percent.  In  this  case  the  initial 
sideband  in  the  simulation  of  Section  2.2  lies  well  within  the  bandwidth  of 
the  cavity.  The  loss  at  the  sideband  frequency  is  2.5  times  that  of  the 
main  line.  It  is  found  that  these  losses  are  not  sufficient  to  suppress 
all  sideband  activity  and  restore  full  extraction  efficiency.  Evidently  a 
cavity  with  strong  dispersion  is  required  for  the  PEL  parameters  studied 
here. 

In  the  case  of  the  4  percent  bandpass  filter,  the  operating 
wavelength  stays  fixed  at  the  startup  wavelength.  The  line  does  not  chirp 
due  to  the  rather  steep  rise  in  loss  to  either  side  of  the  main  line.  In 
the  case  of  the  edge  filter,  the  losses  at  first  rise  more  slowly  with 
increasing  wavelength,  allowing  the  line  to  chirp  about  0.6  percent  in 
wavelength.  This  operating  point  is  more  nearly  the  point  of  peak 
saturated  gain.  In  fact,  in  thiB  case  the  laser  picks  am  operating  point 
with  about  13  percent  output  coupling  rather  than  the  nominal  10  percent. 

The  best  choice  of  a  filter  shape  depends  on  several  factors. 

Modeling  predicts  that  many  tapered-wiggler  PEL  designs  will  require 
substantial  wavelength  chirp  during  the  transition  from  small-signal  to 
saturated  operation,  so  in  that  case  the  filter  function  must  be  configured 
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Figure  2-9.  Introduction  of  Wavelength  Selectivity  Eliminates  Loss  of 
Extraction  Efficiency  Due  to  Sidebands. 


2-25 


to  allow  for  the  chirp.  The  wiggler  taper  used  in  this  example  is 
especially  designed  to  minimize  chirp  requirements;  it  will  start  up 
without  chirp  at  an  output  coupling  of  10  percent,  although  a  small  amount 
of  chirp  to  longer  wavelengths  is  required  for  startup  with  larger 
outcoupling.  Another  consideration  is  loss  induced  by  the  chirp.  In  the 
edge  filter  example  shown  in  Figure  2-8,  the  laser  has  considerably  higher 
loss  after  the  chirp.  This  is  desirable  only  if  the  additional  losses  are 
recoverable  as  outcoupled  power. 
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Section  3 


NAVEUNGTH-SEIMCTIVITY  OPTIONS 

Wavelength-dispersive  optical  elements  may  be  required  in  the  laser 
cavity  for  suppression  of  the  sideband  instability.  This  section  discusses 
several  schemes  which  provide  adjustable  selectivity  in  the  desired 
wavelength  range  of  about  ±1  percent,  i.e.,  in  a  range  of  around  5  to  10 
nm.  The  wavelength-selective  element  must  provide  sufficient  losses  to 
suppress  the  sideband  evolution  without  adding  losses  to  the  cavity  which 
prevent  startup  at  the  design  wavelength.  The  selective  optics  must  add 
very  little  wave  front  distortion,  around  X/100  per  surface.  This  implies 
highly  homogeneous  materials,  especially  for  transmissive  elements,  and 
very  high  quality  optical  surfaces. 

The  approaches  considered  include  thin-film  dielectric  coatings  on 
the  end  mirrors,  birefringent  filters,  and  gratings,  as  well  as  a  novel 
approach  based  on  optical  mode  size  variation  with  wavelength  using  a 
dispersive  lens  in  the  cavity.  We  have  considered  designs  which  would  be 
suitable  for  near-term  visible-wavelength  PEL  oscillator  experiments,  but 
note  that  the  most  promising  candidates  for  high-power  systems,  which 
preclude  use  of  transmissive  optics  within  the  cavity  due  to  cooling 
requirements,  are  dielectric  mirror  coatings  and  gratings.  While  not 
suitable  for  high-power  designs,  the  birefringent  filter  has  the  property 
to  also  serve  as  the  cavity  outcoupler.  Thus  any  additional  cavity  losses 
induced  by  the  filter  following  frequency  chirp  of  the  PEL  cam  be  included 
in  the  outcoupled  power,  am  importamt  efficiency  consideration. 

3.1  THIN-FILM  MIRROR  COATINGS 

Thin-film  dielectric  coatings  are  often  used  on  cavity  end  mirrors  to 
achieve  high  reflectivity  at  the  design  frequency.  Such  mirrors  generally 
are  highly  reflective  only  within  a  fairly  narrow  band  of  wavelengths  and 
cam  therefore  be  used  for  wavelength  selectivity.  The  reflectivity 


_  • 
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profiles  of  various  realizable  dielectric  stack  designs  are  calculated  and 
evaluated  for  suitability  as  wavelength-selective  elements  in  an  FEL. 

A  computer  algorithm  was  written  to  calculate  the  ideal  reflectance 
properties  of  thin-film  periodic  multilayer  stacks/ 3-1  ^  Each  period  of 
the  stack  consisted  of  either  two  or  four  materials  with  alternating  high 
and  low  refractive  indexes.  The  refractive  index  was  assumed  homogeneous 
and  the  coating  thickness  was  assumed  uniform  across  the  aperture  of  the 
mirror. 

It  should  be  noted  that  the  model  does  not  account  for  the  following 
real  aspects  of  multilayer  stacks/3-2  to  3-4 )  The  polycrystalline 
structure  of  the  high  index  material  will  lead  to  light  scattering  and 
absorption  of  light.  Hiis  effect  may  be  modeled  by  giving  its  refractive 
index  an  imaginary  component  corresponding  to  the  absorption  at  a  given 
wavelength.  Due  to  the  paucity  of  experimental  data  for  the  coatings  of 
interest,  this  effect  was  ignored.  For  realistic  cases,  the  shape  of  the 
functions  used  in  the  sideband  calculations  of  Section  2.3  will  not  change 
significantly,  and  the  conclusions  of  that  Section  will  still  hold.  The 
coating  designs  considered  were  purposely  limited  to  physical  thicknesses 
of  *5  /a  or  less  since  errors  due  to  cumulative  film  nonuniformity  across 
the  effective  mirror  aperture  degrade  the  Strehl  ratio  significantly  in 
thick  stacks.  For  example,  a  cumulative  nonuniformity  of  only  2  percent  in 
a  4  Mm  thick  stack  causes  an  effective  figure  error  of  X/s/3-5^  We  note 
that  assuming  perfect  substrates  and  a  coating  figure  error  of  X/28  causes 
the  on-axis  intensity  of  a  mirror  to  decrease  20  percent.  Such 
considerations  dictate  control  of  thickness  uniformity  to  below  the  0.5 
percent  level. 

The  purpose  of  the  calculation  was  to  determine  which  features  of 
thin-film  stacks  were  most  significant  for  wavelength  selectivity  in  the 
FEL.  Various  designs  and  combinations  of  materials  were  used  in  the 
calculations,  The  stack  had  either  one  or  two  high/low  index  of  refraction 
combinations  per  period.  The  high/low  combinations  were  either  Al~03/Si02, 


Zr02/Si02,  or  Tio2/Si02,  with  refractive  index  ratios  of  1.126,  1.499,  and 
1.8,  respectively.  The  bulk  refractive  indexes  were  used  in  all  of  the 
coating  designs,  due  to  lack  of  precise  knowledge  of  the  actual  thin-film 
refractive  index  for  all  but  one  of  the  materials.  The  thin-film 
refractive  index  values  will  always  be  less  than  the  bulk  values,  and  thus 
the  reflectivities  of  an  actual  coating  may  be  expected  to  be  slightly  less 
than  these  theoretical  curves. 

Results  of  varying  the  refractive  index  ratio  for  three-quarter  wave 
layers  are  presented  in  Pigure  3-1.  For  a  fixed  number  of  periods,  N,  the 
maximum  reflectivity  R^,,  increases  with  the  refractive  index  ratio,  n^n^, 
as  does  the  wavelength  range  over  which  the  reflectivity  decreases  by  less 
than  1  percent  (the  high  reflectivity  bandwidth).  The  designs  shown  in  the 
figure  have  l^,T  occurring  at  0.5  tm.  It  should  be  noted  that  for  the 
physical  thickness  used,  the  A1203  design  does  not  provide  high 
reflectivity  and  may  be  useful  only  if  it  also  serves  as  the  outcoupler. 

Figure  3-2  demonstrates  that  increasing  N  for  quarter-wave  designs 
increases  both  R[nl.r  and  the  bandwidth.  It  may  be  observed  from  the  results 
of  Section  2.3  that  simple  quarter-wave  designs  with  a  high  reflectivity 
fail  to  provide  narrow  enough  bandwidths  for  wavelength  discrimination. 

The  designs  using  3X/4  thick  layers  in  each  period  shown  in  Figure  3-1  are 
superior  in  this  respect.  Figure  3-3  illustrates  this  design  feature  for  a 
fixed  number  of  periods  of  the  ZrO2/Si02  combination,  which  is  the  most 
likely  candidate  for  PEL  applications  (because  of  its  damage  resistance 
properties).  Either  of  the  3X/4  or  5X/4  designs  appear  to  provide  the 
desired  sideband  suppression  characteristics .  Since  both  damage 
sensitivity  and  figure  error  increase  with  the  physical  thickness  of  the 
stack,  the  3X/4  coating  would  probably  be  chosen. 

Figure  3-4  demonstrates  the  kind  of  edge  filter  type  coating  designs 
that  may  be  achieved  by  choosing  more  complicated  structures  for  each 
coating  period.  In  this  case  am  alternating  high/low/high/low  combination 
was  employed  with  the  second  high/ low  combination  corresponding  to  a 
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Figure  3-3.  Dependence  on  Reflectivity  on  Film  Thickness  for 
ZrO„  +  SiO_  Ten  Period  Stack. 
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Figure  3-4.  Five  Period  Edge  Filter  Design  for  Two  Layer 


Thicknesses  of  Zr0„  +  SiO„. 
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different  wavelength  of  maximum  reflectivity.  As  can  be  observed  in  the 
figure,  the  edge  filter  slopes  can  be  made  very  sharp  so  such  coating 
designs  should  completely  suppress  sideband  development. 

The  damage  threshold  claims  for  edge  filters  are  just  as  high  as  for 
the  best  "hard"  mirror  coatings  (i.e.,  a  claim/3-6)  of  5  MW/ cm 2  CW  at  1.06 
tun).  Most  coating  designs  include  a  X/2  overcoat  of  Si02,  since  thia  has 
been  shown  experimentally*3-7)  to  increase  the  damage  threshold  by  at  least 
a  factor  of  2.  The  effect  of  this  additional  coat  on  the  maximum 
reflectivity,  bandwidth  and  slope  is  insignificant  as  far  as  PEL 
performance  and  sideband  suppression  is  concerned. 

3.2  BIREPRHKSDfT  FILTERS 

Birefringent  filters,  also  called  Lyot  filters,  sure  based  on  the 
natural  birefringence  of  certain  materials  such  as  crystalline  quartz , 
which  is  the  only  material  considered  here.  Frequency  tuning  is  provided 
based  upon  the  differing  phase  retardation  of  ordinary  and  extraordinary 
waves  which  pass  through  the  crystal.  A  plane  polarized  input  wave  will 
become  elliptically  polarized  in  a  manner  dependent  on  the  input 
wavelength,  because  of  the  dispersion  of  both  ordinary  and  extraordinary 
refractive  indexes.  Filtering  occurs  because  the  FEL  provides  gain  in  only 
one  plane  (assuming  a  plane-polarized  wiggler),  and  light  which  suffers  a 
90  degree  rotation  in  a  birefringent  plate  tilted  at  near  Brewster's  angle 
will  undergo  considerable  loss  by  reflection  ("17  percent  per  surface). 
Figure  3-5  defines  the  plate  geometry  relative  to  the  FEL.  Light  which  hats 
a  wavelength  such  that  it  is  not  rotated  out  of  the  wiggler  plane  (i.e.,  it 
remains  a  TM  wave)  suffers  virtually  no  loss.  Thus  these  devices  are 
nearly  ideal  filters  in  that  their  transmission  function  has  a  theoretical 
maximum  of  unity.  Their  demonstrated  high  damage  resistance  is  an 
additional  advantage . * 3-8 ) 

The  effective  rotation  of  the  plane  of  polarization  as  a  function  of 
X  is  described  in  terms  of  a  2  x  2  Jones  matrix  which  is  a  function  of  9, 
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Figure  3-5.  Orientation  of  Tilted  Birefringent  Plate  in  FEL  Cavity 
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the  plate  tilt  angle  with  respect  to  the  system  optic  axis,  and  0,  a 
rotation  of  the  optic  axis  within  this  plane  (see  Figure  3-5).  For  low 
loss  6  is  chosen  equal  to  Brewster's  angle,  the  center  wavelength  is 
chosen  by  variation  of  0.  The  general  form  of  the  transmission  function 
for  the  TM  component  of  the  incident  wave  is^ 3-9  ^ 


T(0,0)  =  1  -  A(0,0)  ein 


afOjX)  -  0  (X) 


[3-D 


where  A  is  the  depth  of  modulation  of  the  transmission  function  and  0e  and 
0o  are  the  wavelength -dependent  phase  shifts  inputed  by  the  birefringent 
plate.  When  the  argument  of  the  sin  function  is  a  multiple  of  v,  the 
transmission  function  is  unity.  The  thickness  of  the  birefringent  plate 
and  its  angle  0  determine  the  period  of  T (0,0)  and  its  wavelength  of 
maximum  transmission. 


Figure  3-6  illustrates  the  tunability  of  three  filters  with 
thicknesses  of  400,  800  and  1200  fm.  These  filters  are  inserted 
intracavity  at  Brewster's  angle.  With  0-45  degrees,  the  maximum  in  the 
transmission  function  occurs  at  0.55  ftm.  Other  wavelength  values  are 
easily  obtained  by  simple  rotation.  The  depth  of  modulation  is  seen  to 
depend  on  this  angle  as  well.  Figure  3-6  demonstrates  the  increase  in 
slope  obtained  by  going  to  thicker  plates. 

High  quality  plates  axe  readily  available  in  crystalline  quartz  at 

all  the  design  thicknesses  considered  here.  Typical  scattering  losses 

below  0.1  percent  per  surface  axe  desired  and  can  be  manufactured .  A 

surface  flatness  approaching  X/100  can  be  achieved,  so  introduction  of 

additional  figure  error  is  not  a  problem.  Power  levels  exceeding  500 
2 

MW/cm  in  a  Q  switched  pulse  train  at  a  pulse  repetition  rate  of  10  pps 
have  been  achieved  in  intracavity  Nd  t Glass  laser  systems  and  no  spectral 
shift  due  to  thermal  effects  was  observed . ^ 3-8  ^  it  should  be  noted  that 
the  latter  effect  is  second  order  for  a  birefringent  filter  as  compared  to 
a  solid  Fabry-Perot. 
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Figure  3-6.  Birefringent  Plate  Transmission  Functions  for  Various 
Plate  Thicknesses,  t.  (Dotted  curve  is  offset  for 
Clarity. ) 


A  feature  of  this  filter  is  the  possibility  of  providing  wavelength- 
dependent  outcoupling.  If  one  side  of  the  birefringent  plate  is 
antireflection  coated  and  the  PEL  cavity  is  a  ring,  the  filtered  light  is 
in  a  single  reflected  bean.  Furthermore,  if  the  plate  is  tilted  slightly 
off  Brewster's  angle,  it  can  provide  the  primary  means  of  outcoupling  for 
the  cavity,  with  the  outcoupled  light  directed  in  the  same  beam  as  that 
reflected  by  the  filter.  If  the  PEL  chirps  during  the  transition  from  the 
small-signal  to  saturated  regimes,  as  is  generally  the  case  for  tapered- 
wigglers,  the  saturated  outcoupling  could  exceed  the  outcoupling  during 
startup.  This  is  useful  for  minimizing  the  startup  time. 

3.3  ANGUIARLY— DISPERSIVE  EUSHENTS 

Angularly— dispersive  elements  include  both  gratings  and  prisms.  Such 
elements  can  provide  wavelength  selectivity  in  an  PEL  cavity  by  dispersing 
optical  sidebands  so  that  they  lose  overlap  with  the  gun  medium  and  are 
not  amplified.  In  the  case  of  the  prism,  the  wavelength  discrimination  is 
essentially  fixed  by  the  dispersion  of  the  refractive  index  of  the  prism 
material  and  the  beam  size.  Since  a  prism  would  be  expected  to  handle 
rather  high  power  levels  in  transmission,  fused  quartz  would  likely  be 
required.  Operation  at  a  tilt  angle  near  Brewster's  angle  would  be 
necessary  to  minimize  insertion  losses.  With  these  design  constraints,  a 
frequency  resolution  of  AX/X  «  0.05  percent  would  be  expected  for 
conditions  of  a  visible-wavelength  PEL  such  as  the  design  considered  in 
Section  2.2.  This  may  actually  be  too  severe  for  PEL  requirements. 
Diffraction  gratings  used  in  reflection  offer  more  design  freedom  since  the 
resolution  is  set  by  the  grating  period.  For  these  reasons,  the  following 
design  analysis  focuses  on  gratings. 

For  a  grating  to  operate  as  a  sideband  suppressor,  the  dispersion  of 
the  grating  must  be  large  enough  to  misalign  the  sideband  component  and 
small  enough  that  all  frequency  components  of  the  desired  micropulses  are 
still  aligned  in  the  cavity.  This  is  possible,  since  the  transform- limited 
spectral  width  of  the  micropulse  is  much  less  than  the  sideband  shift. 


Another  requirement  is  efficiency.  If  the  grating  is  not  the  output 
coupler,  the  grating  efficiency  should  be  99  percent  or  greater.  Another 
consideration  Which  may  limit  the  selectivity  is  potential  cavity 
misalignment  due  to  chirp  during  the  startup  phase . 

3.3.1  Grating  Design  for  FEL  Applications 

The  calculations  of  Section  2  indicate  that  about  1  percent 
wavelength  discrimination  eliminates  sideband  development.  As  a  design 
example,  consider  a  grating  used  in  Littrow  or  Hautocollimation"  mounting 
as  shown  in  Figure  3-7,  which  we  will  configure  to  suppress  a  sideband  of 
wavelength  offset  AX/X  of  1.5  percent.  For  this  geometry  the  angle  of 
incidence  equals  the  blaze  angle  and  the  angle  of  diffraction.  The 
Rayleigh  criterion  for  resolution  of  the  grating  is  determined  by  the 
optical  path  difference  ( in  waves  of  light )  between  the  extremities  of  the 
grating*  3-10  > 


2W  sin0  ' 


[3-2] 


where  W  is  the  grating  width  and  0  is  the  angle  of  incidence,  we  assume 
that  the  beam  size  in  the  cavity  requires  a  grating  with  a  width  of  "2  cm. 
The  ideal  blazed  grating  will  have  triangular  grooves,  and  the  angle  that 
the  primary  facet  makes  with  respect  to  the  grating  plane  (the  blaze  angle) 
will  be  chosen  to  be  the  same  as  the  angle  of  incidence  (see  Figure  3-7). 
Each  groove  is  formed  so  that  independently  (by  geometric  optics)  it  acts 
like  a  small  mirror  tilted  to  redirect  the  light  in  the  direction  of  a 
chosen  diffracted  order.  For  the  FEL  described  in  Section  2.2,  the  grating 
will  have  a  blaze  angle  of  ”0.05  degrees.  The  line  spacing  will  be  ”0.3 
mm,  an  extremely  modest  requirement. 


The  groove  profile  determines  the  amount  of  energy  diffracted  into 
any  individual  order.  We  define  the  absolute  efficiency  to  be  the  fraction 
of  the  incident  power  Which  is  distributed  in  a  single  diffracted  order. 
This  is  the  efficiency  quoted  herein  unless  otherwise  noted.  The  relative 
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Figure  3-7.  Illustration  of  Diffraction  Grating  Used  in 
Littrow  Mounting. 


efficiency  is  defined  as  the  fraction  of  the  incident  power  which  would  be 
distributed  in  a  single  diffracted  order  if  the  grating  surface  had  perfect 
reflectivity.  Thus  the  relative  efficiency  is  a  measure  of  the  groove 
efficiency.  The  periodicity  and  aperture  function  of  the  grating  influence 
only  the  resolution  of  the  grating,  and  because  of  the  modest  resolution 
requirements  of  the  PEL,  manufacturing  imperfections  in  these 
characteristics  are  less  important  than  control  of  the  groove  geometry. 
Since  absolute  efficiencies  for  intracavity  application  of  gratings  for  the 
PEL  should  exceed  99  percent,  we  may  immediately  rule  out  symmetric  groove 
profiles  which  can  be  at  best  only  50  percent  efficient  at  near  normal 
incidence  angles. 

3.3.2  Production  of  Concave  Biased  Gratings 

Pigure  error  budget  considerations  for  the  PEL  dictate  a  minimum 
number  of  intracavity  optical  elements.  Although  planar  gratings  would  be 
of  interest  in  ring  cavity  designs,  concave  gratings  are  useful  in  linear 
cavities  since  they  combine  both  focusing  and  dispersion  in  a  single 
element. ^ 3-11 )  Thus  the  design  and  manufacture  of  concave  gratings  will  be 
addressed . 

Historically,  most  gratings  were  manufactured  by  diamond  ruling  of  an 
aluminum  or  gold  coated  planar  glass  substrate.  This  technique  will  not 
yield  gratings  with  the  desired  efficiencies  for  PEL  applications,  since  it 
is  very  difficult  to  control  groove  profiles  for  blaze  angles  less  than  0.5 
degrees  due  to  tool  geometry  and  other  factors. ( 3-10 )  in  addition,  to 
achieve  the  required  reflectivities  at  near  normal  incidence  will  require 
multilayer  dielectric  coated  gratingB,  which  are  difficult  to  fabricate  by 
ruling  techniques.  If  gratings  are  ruled  on  a  spherical  substrate  using 
the  same  engine  as  employed  for  planar  gratings,  the  groove  spacing  and 
facet  angle  will  not  be  uniform  with  respect  to  the  spherical  grating 
surface.  If  conventional  ruling  techniques  must  be  employed,  efficient, 
aberration-free  gratings  with  a  constant  groove  angle  may  be  made  by 
cylindrically  bending  a  diffraction  grating  ruled  on  a  plane 


surf ace . ( 3~ i0 )  However,  due  to  the  inconvenience  of  manufacture  of  ruled 
concave  gratings  and  the  lesser  efficiency  of  such  gratings,  interference 
gratings  can  be  an  attractive  alternative. 

Hie  simplest  technique  of  creating  interference  gratings  consists  of 
coating  a  blank  substrate  with  photoresist,  exposing  it  to  a  fringe  pattern 
and  developing  it.  The  intensity  variation  translates  into  a  solubility 
variation  across  the  photoresist.  After  development  with  a  suitable 
solvent,  the  remaining  photoresist  has  the  proper  spatial  frequency 
variation.  A  reflective  coating  can  then  be  applied.  The  control  of  the 
thickness  and  uniformity  of  the  photoresist  coating  is  of  paramount 
importance  for  the  overall  quality  of  the  resulting  grating.  Direct 
application  of  interference  grating  techniques  at  normal  incidence  produces 
gratings  with  symmetric  groove  profiles. 

A  blazed  grating  can  be  produced  by  a  simple  interference  scheme 
which  involves  tilting  the  grating  blank.  A  convex  mirror  is  placed  behind 
the  grating  blank  and  centered  on  the  desired  focal  point  of  the  grating. 

A  standing  wave  pattern  is  formed  at  the  tilted  grating  blank  surface  which 
produces  fairly  flat  facets  whose  normals  all  point  to  the  same  focal  point 
(the  condition  which  provides  uniform  efficiency  and  constant  blaze 
wavelength  over  the  grating  aperture).  The  wavelength  of  the  laser  used  to 
set  up  the  standing  wave  interference  pattern  will  determine  the  blue 
wavelength  of  the  grating.  This  technique  has  yielded  groove  profiles  of 
sufficient  quality  to  produce  99  percent  relative  efficiency  at  a  blue 
angle  of  10  degrees. ^ 3-12 ) 

For  FEL  applications,  somewhat  shallower  blue  angles  will  likely  be 
required.  The  blaze  angle  of  an  interference  grating  can  be  reduced  by  ion 
beam  etching.  This  technique  hu  produced  high  quality  gratings  with 
angles  near  0.5  degrees  and  hu  demonstrated  blue  angle  reduction  by 
factors  of  3  to  10.^  3-13  ^  Gratings  blued  for  wavelengths  between  0.5  im 
and  5  Mm,  with  efficiencies  of  up  to  90  percent,  have  been  prepared  with 
this  technique . ^ 3_1* ) 


Another  technique  of  interest  for  FEL  applications  uses  the 
preferential  chemical  etching  properties  of  crystalline  silicon  to  produce 
a  blazed  grating  following  exposure  of  a  photoresist  surface  layer  to  a 
symmetrical  interference  pattern.  A  suitable  photoresist  is  deposited  on  a 
Si  crystal  oriented  with  its  <110>  plane  at  an  appropriate  angle  to  the 
surface.  The  <no>  plane  etches  at  a  greater  rate  than  the  <lll>  plane 
resulting  in  a  triangular  groove  profile.  Thus  the  blaze  angle  is 
determined  by  the  angle  of  the  crystal line  axis  relative  to  the  surface, 
rather  than  by  the  photoresist  exposure  process.  The  resulting  surface  may 
be  coated  with  a  dielectric  stack  for  good  efficiency.  Features  to  10-4 
cm  can  be  created  in  this  manner.  Typical  rms  surface  roughness  is  2  to  5 
Angstroms.  The  smoothness  of  the  surface  provides  low  scattering  loss  and 
high  damage  resistence,  two  attributes  required  of  FEL  optics. 

3.3.3  Practical  Performance  of  Interference  Gratings 

The  best  interference  and  ruled  gratings  produce  diffracted  wave 
fronts  of  high  optical  quality,  but  interference  gratings  suffer  less  stray 
light  and  scattering  problems.  Scattering  losses  for  interference  gratings 
are  predicted  to  be  less  than  0.2  percent. 

The  theoretical  efficiency  of  a  grating  with  two  dielectric  overcoats 
is  "95  percent.  An  experimental  grating  with  0.015  tm  groove  depth 
(approximately  the  value  required  for  the  FEL)  has  been  made  with  87 
percent  efficiency;  the  loss  in  efficiency  was  believed  to  be  due  to  the 
inability  to  control  the  blaze  angle  over  the  entire  grating 
aperture/3-15^  Ho  re  recently,  a  grating  with  nearly  99  percent  relative 
efficiency  was  fabricated  using  preferential  etching  in  Si  and  ion 
etching/3-12)  Presumably,  with  suitable  dielectric  coatings  this  type  of 
grating  would  be  suitable  for  FEL  application,  since  most  of  the 
interference  gratings  with  low  blaze  angles  are  first  generation 
prototypes,  one  may  expect  improved  efficiencies  as  experience  in  their 
manufacture  is  gained. 


Highly  damage  resistent  gratings^  3-16 )  are  usually  ruled  in 

substrates  of  pyrex  or  copper  (for  infrared  applications ) .  Typical  CW 

damage  levels  at  1.06  (im  wavelength  are  10  XW/cm2  for  10  sec  (glass 

2  2 

substrate),  while  pulsed  levels  are  140  MW/cm  at  25  nsec,  70  MW/cm  at 

2 

100  nsec,  and  7  Joules/ cm  at  100  nsec  (on  an  electroless  nickel  copper 
substrate).  Thus  circulating  powers  could  be  as  high  as  35  MW  for  100  nsec 
pulses  at  low  repetition  rates.  Thus  CW  damage  would  limit  the  design  of 
Section  2.2  (but  at  1  fm)  to  several  kW  average  circulating  power.  The 
coating  material  for  the  above  figures  is  aluminum.  Overcoating  with  gold 
or  silver  improves  the  reflectivity  but  lowers  the  damage  threshold  by 
2.  Damage  measurements  of  dielectric  coated  gratings  are  not  available, 
but  would  be  expected  to  be  higher,  as  is  the  case  for  mirrors. 

3.4  FREQUENCY  SELECTIVITY  THROUGH  A  DISPERSIVE  OPTICAL  MODE  SIZE 

The  PEL  optical  gain  is  sensitive  to  spatial  characteristics  of  the 
optical  mode  due  to  variation  of  the  electron-photon  overlap  with  mode 
size.  This  property  can  be  exploited  to  provide  a  frequency-dependent 
round  trip  loss  for  the  purpose  of  suppressing  the  sideband  instability. 

The  basic  concept  is  to  use  a  dispersive  lens  in  the  cavity  to  provide  a 
wavelength-dependent  mode  size,  and  to  configure  the  cavity  so  that  the 
available  dispersion  provides  the  degree  of  selectivity  desired. 

A  convenient  property  of  the  PEL  so  far  as  this  goal  is  concerned  is 
the  need  for  near-concentric  cavities.  This  geometry  provides  the  small 
spot  size  at  the  wiggler  and  the  large  spots  at  the  end  mirrors  necessary 
to  avoid  mirror  damage.  The  mode  size  of  the  near-concentric  cavity  is 
extremely  dependent  on  the  mirror  focal  length  (assuming  fixed  cavity 
length),  so  that  only  small  focal  length  changes  are  needed  from  the 
dispersive  element.  The  Rayleigh  range,  Zg,  for  the  focal  region  near  the 
center  of  a  symmetric  near-concentric  cavity  is  a  simple  function  of  the 
distance  between  the  centers  of  curvature  of  the  end  mirrors,  d,  and  is 
given  by 


where  L  is  the  total  cavity  length.  The  Rayleigh  range  determines  the  mode 
size,  and  the  mode  size  determines  the  electron  photon  overlap  and,  hence, 
the  FEL  gain.  An  example  calculation  of  gain  as  a  function  of  distance 
between  the  centers  of  curvature  of  the  end  mirrors  is  given  in  Figure  3- 
8.  Key  parameters  of  this  calculation  are  listed  in  Table  3-1.  It  is 
worth  noting  that  the  overall  shape  of  this  curve  is  dependent  primarily  on 
the  emittance  value  chosen.  Aperturing  of  the  beam  by  the  wiggler  magnets 
is  unimportant  over  the  range  shown. 

The  geometry.  Which  can  provide  variable  selectivity  at  the  1  percent 
wavelength  resolution  of  interest,  is  shown  in  Figure  3-9.  The  combined 
focal  length  of  the  lens  and  mirror  Mg  becomes  more  dispersive  as  their 
focal  lengths  decrease,  while  the  power  of  the  combination  can  remain  fixed 
if  the  spacing  between  them  is  adjusted  properly.  The  mirror-lens 
combination  was  analyzed  with  standard  ABCD  matrix  ray-trace  formalism  to 
determine  the  sensitivity  of  the  combination  to  dispersion.  As  am  example, 
consider  a  combination  which  has  a  positive  fused  silica  lens  with  aim 
focal  length  separated  from  mirror  Mg  by  10  cm.  The  corresponding  convex 
mirror  to  produce  an  effective  focal  length  of  15  m  has  a  focal  length  of 
-0.457  m.  The  sensitivity  of  this  example  to  wavelength  is  illustrated  by 
a  28  cm  change  in  focal  length  for  a  5  nm  wavelength  change .  The  focal 
length  variation  can  be  combined  with  Figure  3-8  for  a  rough  estimate  of 
the  net  gain  as  a  function  of  wavelength,  and  is  given  in  Figure  3-10.  The 
changes  in  effective  focal  length  for  5  nm  excursions  of  the  laser 
wavelength  are  the  correct  size  to  suppress  the  sideband  instability,  and 
still  provide  acceptable  gain  over  the  range  of  wavelengths  desired  for  FEL 
startup .  The  selectivity  of  the  lens  mirror  combination  can  be  adjusted  by 
choosing  different  combinations  of  strengths  of  the  elements. 

The  main  detriment  to  the  use  of  the  lens-mirror  combination  in  the 
FEL  cavity  is  the  effect  on  the  figure  error  budget.  The  introduction  of 
the  lens  into  che  cavity  adds  a  figure  error  Which  is  sensed  twice  each 


Table  3-1 


PARAMETERS  OP  GAIN  CALCULATION  WITH  VARIABLE  NODE  SIZE 


Normalized  Quittance 

0.01  cm  rad 

Higgler  Pull  Gap 

0.45  cm 

Higgler  Length 

5  m 

Higgler  Havelength 

2.02  cm 

Optical  Havelength 

0.5  Mm 

Beam  Energy 

120  MeV 

e  B  X 

max  w 

aw  '  2/3  2 

2  '  mac 

1.33 

Cavity  Length 

60  meters 
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Near-Concentric  Cavity  Geometry  for  Use  of 
Dispersive  Lens  for  Mode  Size  Variation. 
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WAVELENGTH  CHANGE  (nm) 


Figure  3-10. 


Calculation  of  Net  FEL  Gain  Variation  Due  to  Change 
in  Mode  Size. 


round  trip.  The  lens  figure  error  introduced  by  the  two  passes  each  round 
trip  is  correlated  since  a  ray  passes  through  the  same  spot  on  the  lens  on 
each  pus.  For  example,  a  lens  made  to  a  figure  error  of  X/50  introduces 
an  effective  round  trip  figure  error  of  X/25,  which  in  some  cases  may  be 
an  unacceptable  addition  to  the  figure  error  budget. 
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Abstract 


Electron-beam  energy  spectral  measurements  were  made  on  a  tapered-wiggler  free-electron 
laser  amplifier.  A  10  MeV  electron  beam  from  a  traveling-wave  linear  accelerator  inter¬ 
acted  in  a  tapered-wiggler  with  an  intense  10.6um  CO2  laser  beam.  The  electron  spectra 
show  a  4  percent  net  energy  loss  and  a  9  percent  peak  loss.  Measurements  of  electron 
energy  spectra,  extraction  efficiency  as  a  function  of  electron-beam  energy,  and  extraction 
efficiency  as  a  function  of  optical  power  are  presented  and  are  consistent  with  theoreti¬ 
cally  predicted  performance. 


Introduction 


The  experimental  results  reported  here  help  validate  the  concept  of  the  taper-wiggler 
free-electron  laser  (FEL) .  The  device  is  a  candidate  for  a  high-efficiency  tunable  source 
of  coherent  radiation.  The  first  FEL,  demonstrated  at  Stanford,1  had  an  untapered  wiggler. 
The  tapered  wiggler2  differs  from  the  Stanford  wiggler  in  that  the  resonant  electron  energy 
of  the  wiggler  magnet  varies  along  its  length  to  maintain  a  resonant  electron-photon  inter¬ 
action  as  the  electrons  decelerate.  The  resonant  energy  of  the  wiggler  can  be  tapered  by 
varying  the  wavelength  or  amplitude  of  the  periodic  magnetic  field  as  a  function  of  axial 
position.  Electrons  trapped  in  the  ponderomotive  potential  well  formed  by  the  electric 
field  of  the  light  and  the  magnetic  field  of  the  wiggler  decelerate  in  accord  with  the 
resonant  energy  change,  or  taper,  of  the  wiggler.  Adding  taper  to  the  FEL  wiggler  can  pro¬ 
vide  increased  electron-beam  energy  extraction  and  increased  overall  efficiency  at  the 
expense  of  reduced  gain. 

Not  all  of  the  electrons  entering  the  wiggler  are  trapped  in  the  ponderomotive  potential 
and  the  nominal  trapping  fraction  for  this  experiment  is  about  50  percent.  Under  these 
conditions  the  net  deceleration  is  about  half  the  peak  deceleration.  The  corresponding 
energy  spectrum  has  two  peaks  of  roughly  equal  current,  one  near  the  entrance  energy  and 
the  other  at  the  energy  resonant  with  the  wiggler  exit.  The  wiggler  used  in  this  experi¬ 
ment  has  a  9  percent  energy  taper  and  a  nominal  net  electron  energy  extraction  of  about 
4  percent. 

The  highest  net  energy  extraction  previously  reported  by  our  group  is  about  2.5  percent 
indicating  strong  interaction,  but  less  than  optimal  trapping  over  the  length  of  the  9  per¬ 
cent  taper  wiggler  in  our  experiment.3  This  result  is  limited  by  a  combination  of  non¬ 
optimum  overlap  of  the  electron  and  laser  beams,  the  electron  beam  emittance,  and  the  laser 
power.  Detrapping  of  electrons  along  the  wiggler  length  is  evident  in  the  energy  spectra. 
Here  we  report  results  of  electron-beam  extraction  measurements  in  which  the  electron  beam 
emittance  is  reduced  from  our  previous  experiments.  He  again  use  the  9  percent  tapered 
wiggler.  The  electron  beam  size  is  not  limiting,  and  4  percent  net  extraction  is  observed. 
The  results  are  consistent  with  theoretical  predictions  and  give  added  verification  to  the 
tapered  wiggler  concept.  They  are  also  consistent  with  the  results  of  the  Los  Alamos 
tapered  wiggler  experiment  in  which  an  energy  taper  of  7  percent  led  to  about  4  percent  net 
extraction. 4 


Experiment 

Electron  beam  energy  extraction  is  measured  in  the  FEL  amplifier  using  the  configuration 
shown  in  Figure  1.  The  wiggler  is  2.3m  long,  has  97  periods,  and  is  constructed  of  SmCos 
permanent  magnets.  It  has  9  percent  resonant  energy  taper  at  constant  synchronous  phase, 
achieved  by  decreasing  the  period  13  percent  and  the  peak  field  strength  8  percent  along 
the  length  of  the  wiggler  with  fixed  gap.  In  the  experiments  reported  here,  the  CO2  laser 
generates  a  40  ns  optical  pulse,  and  the  linear  accelerator  typically  generates  a  0.5  usee 


Figure  1.  Block  diagram  of  amplifier  experiment. 


electron  beam  macropulae  consisting  of  20  psec  micropulses  generated  at  2.8  GHz.  The  two 
beams  are  timed  so  that  the  shorter  optical  pulse  falls  upon  the  longer  electron-beam  pulse 
in  the  wiggler.  During  the  period  of  overlap,  all  of  the  electron  micropulses  are  subject 
to  the  FEL  interaction.  The  spectrograph  can  follow  the  envelope  of  the  temporal  evolu¬ 
tion  of  the  electron-beam  energy  spectra  on  a  nanosecond  timescale,  but  cannot  resolve 
individual  micropulses.  Most  of  the  optical  beam  cannot  interact  with  the  electrons 
because  the  linac  duty  cycle  is  0.05;  therefore,  the  average  gain  is  low  and  is  not  mea¬ 
sured.  The  RF  linac,  electron-beam  spectrograph,  wiggler  magnet,  and  CO2  laser  used  in  the 
extraction  experiments  have  previously  been  reported  upon  in  detail.3 

The  Boeing  linac  is  a  traveling-wave  radio-frequency  linear  accelerator  capable  of 
acceleration  to  30  MeV.  The  primary  power  is  supplied  by  a  single  20  MK  peak,  2  kw  aver¬ 
age,  S-band  klystron.  A  gridded  electron  gun  provides  macropulse  widths  over  the  range  of 
0.005  to  10  psec.  Since  our  last  reported  results,3  the  gun  was  modified  by  the  addition 
of  a  cowling  intended  to  intercept  electrons  from  the  perimeter  of  the  cathode  and  reduce 
the  linac' s  emittance.  Two  stages  of  bunching  at  the  fundamental  frequency  compress  the 
charge  into  15  to  20  degrees  of  phase  of  the  accelerating  wave  to  achieve  an  energy  spread 
of  less  than  2  percent.  The  output  pulse  train  consists  of  pulses  which  are  approximately 
20  psec  long  and  are  separated  by  the  350  psec  RF  cycle  time.  Peak  micropulse  currents  are 
from.  2  to  5  amps  leaving  the  accelerator  structure,  and  as  described  later,  are  lower  at 
the  wiggler  due  to  losses  in  filtering  and  other  losses  in  transport.  The  full  width 
energy  spread  is  2  percent,  and  the  normalized  emittance  for  100  percent  of  the  charge  at 
the  wiggler,  defined  at  yrxx',  is  about  0.01  cm-rad  in  each  plane.  The  beam  is  spatially 
filtered  to  achieve  this  emittance  and  energy  spread,  and  the  peak  micropulse  currents  are 
typically  50  to  200  mA  at  the  wiggler. 

The  electron-beam  transport  and  optics  system  is  shown  in  Figure  2.  The  system  provides 
achromatic  transport  to  the  FEL  with  adjustable  energy  and  emittance  filtering.  The  beam 


Figure  2.  Electron  beam  transport  system. 

is  switched  into  the  FEL  experiment  line  by  a  dog-leg  translation  system  similar  to  a 
design  suggested  by  Swenson.3  Mechanical  slits  in  the  dog-leg  provide  energy  and  emittance 
selection.  Five  quadrupole  magnets  just  upstream  of  the  wiggler  are  used  to  adjust  the 
focal  position  and  convergence  angle  in  two  planes  of  the  input  beam.  Five  fluorescent 
screens  are  deployed  within  the  wiggler  and  six  screens  are  deployed  in  the  transport  legs 
to  assist  in  pointing  and  focusing  the  electron-beam. 

Two  beamline  modifications  have  been  made  since  earlier  experiments.  They  are  the 
upgrade  of  the  emittance  slits  in  the  dog-leg  and  the  addition  of  an  emittance-measuring 
wire-scanner  just  upstream  of  the  wiggler.  The  new  emittance  filter  consists  of  four  inde¬ 
pendently  movable  jaws  and  replaced  a  fixed  aperture  filter.  The  new  jaws  are  typically 


used  with  gaps  3  times  smaller  than  previously  available.  The  emittance  filter  also  has 
been  moved  about  20  cm  down  the  beamline,  placing  it  exactly  at  the  predicted  symmetry 
point  of  the  dog-leg  where  no  encoding  of  energy  on  position  should  exist.  In  order  to 
move  the  emittance  filter,  the  quadrupole  magnet  previously  located  at  the  symmetry  point 
has  been  replaced  by  a  pair  of  quadrupoles  placed  symmetrically  about  the  emittance  filter. 

A  new  emittance-measuring  wire-scanning  profilometer  is  located  directly  upstream  or 
the  wiggler  and  consists  of  a  wire-shadow  scanner,  a  turnout  magnet,  and  a  stopping  block. 
Transport  from  the  scanner  to  the  stopping  block  is  nearly  100  percent,  and  transport  from 
the  scanner  location  to  the  spectrograph  via  the  wiggler,  when  the  turnout  magnet  is  off, 
also  is  near  100  percent.  This  insures  that  the  emittance  and  spectrograph  measurements 
both  include  the  entire  beam.  Emittance  is  deduced  by  measuring  beam  size  as  a  function 
of  the  strength  of  an  upstream  quadrupole  magnet.®*7  These  measurements  give  the  minimum 
spot  size  and  divergence  angle  of  the  beam,  allowing  the  emittance  to  be  computed. 

The  electron  spectrograph  is  a  12-inch  round  pole  Browne -Beuchner 8  design.  The  focal 
plane  of  the  spectrograph  has  segmented  stopping  blocks  cabled  to  oscilloscope  channels. 
These  allow  temporal  bandwidth  of  approximately  200  MHz  and,  as  used,  have  a  minimum  of 
1  percent  energy  spread  per  oscilloscope  channel.  The  use  of  discrete  stopping  blocks 
limits  the  spectrograph  in  energy  resolution,  but  has  the  advantage  of  allowing  temporal 
resolution  and  relatively  simple  calibration. 

The  CO2  laser  consists  of  an  oscillator-preamplifier-amplifier  chain.  A  low  pressure 
gain  cell  forces  the  oscillator  to  operate  on  a  single  longitudinal  mode.  An  electro- 
optical  switch  is  used  to  slice  out  a  fast  rising  40  ns  pulse  from  the  150  ns  oscillator 
output,  as  shown  in  Figure  3.  The  peak  optical  power  delivered  into  the  wiggler  was 
typically  about  0.5  GW  based  on  independent  measurements  of  the  pulse  shape  and  the  inte¬ 
grated  energy. 
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Figure  3.  COj  laser  chain  and  output  pulseform. 
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Electron  energy  spectrum  measurements  and  modeling 

Electron  energy  gain  or  loss  resulting  from  the  PEL  interaction  is  measured  with  the 
electron  spectrograph.  In  these  experiments,  the  interaction  is  easily  identified  in  the 
time-resolved  spectra  because  the  optical  pulse  has  different  temporal  behavior  than  the 
electron  maeropulse.  Representative  histories  of  the  current  into  the  spectrograph  channels 
during  the  CO2  laser  pulse  are  shown  in  Figure  4.  Immediately  prior  to  the  time  of  laser 
overlap,  most  of  the  current  was  entering  the  spectrograph  channel  labeled  -0.5  percent. 
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Figure  4.  Spectrograph  histories  during  FEL  interaction. 

It  can  be  seen  from  the  traces  that  the  mean  energy  is  varying  on  a  timescale  much  longer 
than  the  laser  pulse,  and  this  variation  can  be  ignored.  During  peak  laser  power,  nearly 
90  percent  of  the  electrons  in  the  nominal  input  channel  are  displaced  to  higher  or  lower 
energies.  A  small  amount  of  current  is  accelerated  to  higher  energy  in  the  +1.5  percent 
channel,  and  the  greater  fraction  is  decelerated  into  the  lower  energy  channels,  down  to 
-9  percent.  The  temporal  histories  of  the  channels  differ  greatly  reflecting  the  nonlinear 
nature  of  the  electron  trapping  as  a  function  of  optical  power. 

An  electron  beam  spectrum  measured  at  the  time  of  peak  CO2  intensity  is  shown  in 
Figure  5.  There  is  a  4  percent  shift  in  the  average  electron  energy  assuming  the  current 
is  evenly  distributed  within  each  spectrograph  channel.  The  maximum  extraction  for  any 
electron  is  9  percent.  Superimposed  upon  that  data  is  the  electron  beam  spectrum  prior  to 
injection  of  the  CO2  beam.  The  curves  can  be  considered  as  output  and  input  spectrum, 
respectively,  because  the  macropulBe  current  and  spectral  content  do  not  change  signifi¬ 
cantly  on  several  nanosecond  timescales.  The  photon  pulse,  not  shown,  has  a  peak  power  of 
about  0.5  GW  and  -the  effective  power  could  be  less  because  imperfect  optical  beam  quality 
(Strehl  ratio  <1)  can  only  degrade  the  interaction. 

The  electron  spectra  are  taken  with  the  spectrograph  channels  connected  in  pairs  with 
combined  1  percent  energy  acceptance  over  the  range  of  18.2  to  19.0  MeV,  and  2  percent 
acceptance  elsewhere.  For  any  channel,  the  uncertainty  in  current  at  the  time  of  maximum 
photon  flux  is  less  than  25  percent.  An  indication  of  the  reliability  of  the  measurements 
is  that  the  sum  of  the  measured  channel  currents,  which  should  not  vary,  typically  differ 
"V  *‘»«r  than  10  percent  between  interaction  and  non-interaction  traces.  For  the  data  shown 
in  Figure  5,  the  peak  micropulse  current  at  the  time  of  the  interaction  is  about  160  mA. 

Using  the  input  spectrum  of  Figure  5,  an  output  electron  energy  spectrum,  including  the 
interaction,  has  been  calculated  by  direct  integration  of  the  equations  for  electron  energy 
loss  and  electron  phase  in  the  ponderomotive  potential  of  the  FEL  interaction. 2  In  the 
calculation,  it  is  assumed  that  the  electron-beam  and  500  MW  optical  beams  are  optimally 
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Figure  5.  Measured  input  and  output  electron  beam  spectra  for  FEL  interaction. 


focused  and  coaligned,  and  the  optical  beam  is  diffraction-limited.  A  normalized  emittance 
of  0.023  cm-rad  at  y  »  37  is  assumed  for  both  planes.  The  effect  of  emittance  is  included 
in  the  calculation  in  a  two-step  process.  First  the  optical  electric  (E)  fields  and  wig- 
gler  matnetic  (B)  fields  experienced  by  each  electron  are  computed  as  a  function  of  axial 
position,  including  the  off-axis  motion  but  ignoring  the  small  effect  of  the  FEL  on  the 
trajectory.  Then  the  energy  loss  is  determined  in  a  one-dimensional  integration  of  the 
energy  and  phase  equation  for  electrons  in  the  ponderomotive  potential  well  of  the  FEL. 

This  is  done  using  the  previously  computed  E  and  B  fields  experienced  by  each  electron. 

The  resulting  theoretical  electron  spectrum,  shown  as  the  dashed  line  in  Figure  5,  corres¬ 
ponds  to  a  net  extraction  of  4.0  percent  and  is  in  excellent  agreement  with  the  data.  A 
net  extraction  of  4.7  percent  is  predicted  for  the  same  parameters,  except  with  zero  emit¬ 
tance.  This  shows  that  for  these  parameters  the  extraction  is  only  weakly  dependent  on  the 
emittance.  The  emittance  used  to  give  a  theoretical  curve  matching  the  data  was  about 
twice  the  measured  value.  This  choice  of  emittance  may  roughly  compensate  for  simplifying 
assumptions  used  in  the  theoretical  model  which  would  otherwise  lead  to  an  overestimate  of 
net  extraction.  These  assumptions  are  that  the  CO2  laser  is  diffraction-limited  and  that 
the  electron  and  optical  beams  are  optimally  pointed  and  focused. 

Extraction  as  a  function  of  input  energy  has  been  measured  over  a  range  of  input  ener¬ 
gies  from  18.0  to  19.2  MeV.  This  data  is  shown  in  Figure  6.  The  experimental  data  is 
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Figure  6.  Extraction  as  a  function  of  input  energy. 


upshifted  1.5  percent  in  energy  to  allow  for  a  small  uncertainty  in  spectrograph  calibra¬ 
tion.  Peak  extraction  is  observed  near  19  MeV  for  our  wiggler.  When  the  electron  energy  is 
detuned  below  18.4  MeV,  net  acceleration  of  the  electrons  is  seen.  Electron  deceleration 
and  energy  extraction  is  observed  over  a  range  of  3  percent  in  energy.  This  implies  that 
net  gain  is  present  over  an  optical  bandwidth  of  twice  the  energy  bandwidth  or  about  6  per¬ 
cent.  The  data  is  in  agreement  with  the  theoretically  predicted  curves  generated  assuming 
normalized  emittances  of  0  and  0.023  cm-rad,  perfect  focusing,  perfect  alignment,  and  a 
diffraction-limited  optical  beam.  The  points  shown  do  not  include  all  the  data  taken. 

With  misalignment,  poor  focusing,  or  other  problems,  it  is  always  possible  to  achieve 
results  in  which  the  magnitude  of  the  extraction  is  too  low,  but  it  is  not  possible  to 
achieve  extraction  results  that  are  artifically  high  except  by  actual  sieasurement  errors. 

Electron  energy  extraction  has  been  measured  as  a  function  of  laser  power  and  is  shown 
in  Figure  7.  Extraction  is  seen  tc  increase  with  laser  power.  The  data  is  taken  from 
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Figure  7.  Extraction  as  a  function  of  laser  power  in  the  wiggler. 


three  time-resolved  electron  spectral  records.  Extraction  is  found  as  a  function  of  time 
and  using  laser  pulse  histories,  the  extraction  as  a  function  of  laser  power  is  deduced. 

The  data  is  in  general  agreement  with  theoretical  predictions.  Again,  perfect  focusing  and 
alignment  are  assumed  for  the  theoretical  prediction,  and  a  finite  0.023  cm-rad  normalized 
emittance  is  used  as  a  rough  estimate  of  all  the  effects  of  nonideal  preparation  of  both 
beams.  When  the  extraction  is  predicted  as  a  function  of  laser  power  for  zero  emittance, 
the  theoretical  results  parallel  those  shown  in  Figure  7,  but  the  extraction  is  about 
1  percent  larger. 

Summary 

.*.r.  exoeriment  verifying  the  performance  of  a  tapered-wiggler  FEL  was  described  in  which 
electron-beam  extraction  matching  theoretical  predictions  was  observed.  Deceleration  of 
individual  electrons  by  an  amount  approximately  equal  to  the  energy  taper  of  the  wiggler  is 
observed  as  in  our  earlier  experiments,  a  net  deceleration  of  4  percent  is  observed,  and 
agreement  is  found  between  predicted  and  observed  energy  spectra.  Electron-beam  extraction 
efficiency  is  measured  as  a  function  of  electron  beam  energy  and  CO?  laser  beam  power  and 
the  results  are  consistent  with  predictions  of  performance  at  the  theoretical  limitations. 
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APPENDIX  B 


Electron-Beam  Quality  Requirements  for 
Tapered -Wiggler  Pree-Electron  Lasers 

D.C.  Quimby  and  J.M.  Slater 

Mathematical  Sciences  Northwest,  Inc. 

2755  Northup  Way,  Bellevue,  Washington  98004 

Abstract 

A  general  approach  to  optimization  of  the  free-electron  laser 
interaction  is  used  to  develop  the  scaling  of  emittance  and  energy  spread 
requirements  for  tapered  wigglers  optimized  for  highest  optical  gain  at 
fixed  e-beam  energy  extraction.  The  required  e-beam  properties  for  a  high 
extraction  oscillator  are  found  to  be  quite  stringent  at  visible 
wavelengths,  but  state-of-the-art  accelerators  should  be  sufficient.  The 
applicability  of  various  methods  of  emittance  acceptance  enhancement  is 
examined.  One  very  promising  option  is  a  magnet  canting  scheme  for 
providing  two-plane  focusing  in  planar  wigglers.  Two-plane  focusing 
relaxes  the  severe  emittance  requirement  resulting  from  the  need  to 
maintain  spatial  overlap  between  the  optical  beam  and  the  free-expanding 
e-beam.  In  addition,  options  for  adjusting  various  system  parameters  for 
enhanced  emittance  acceptance,  at  reduced  gain  per  unit  current,  are 
explored . 


I •  Introduction 


The  tapered-wiggler  free-electron  laser  (PEL)  concept  nay  lead  to 
development  of  a  high  efficiency,  tunable  laser.  The  large  e-beaa  energy 
extraction  possible  with  wiggler  tapering  may  lead  to  high  efficiency 
systems  through  reduction  of  e-beam  energy  recovery  or  recirculation 
requirements.  A  key  technology  issue  in  applying  this  technology  at  short 
photon  wavelengths  concerns  the  capability  of  linear  accelerators  to 
produce  electron  beams  with  the  small  emittance  and  energy  spread  required 
by  the  PEL.  The  limited  e-beam  power  available  with  small  emittance  and 
energy  spread,  together  with  the  high  optical  power  required  for  energy 
extraction,  limit  the  tapered-wiggler  single-pass  gain.  The  motivation  of 
this  paper  is  to  gain  a  proper  understanding  of  the  optimization  and 
scaling  of  the  tapered-wiggler  PEL  gain  including  limitations  due  to 
emittance  and  energy  spread,  to  aid  in  development  of  a  high  efficiency, 
high-power  PEL  at  visible  wavelengths.  Emittance  requirements  for 
untapered  wigglers  driven  by  linear  accelerators  have  previously  been 
developed  by  Smith  and  Hadey^1^  and  Dattoli,  et  al . ( 2  ^  In  addition,  Madey 
has  discussed  the  emittance  requirements  for  storage-ring  driven  FELs . ^ 3  ^ 

The  general  approach  to  optimization  of  the  PEL  interaction  strength 
described  in  Section  II  is  used  to  develop  the  scaling  of  emittance  and 
energy  spread  requirements  for  optimized  planar  wigglers.  Emittance  and 
energy  spread  acceptance  limits  are  defined  based  in  part  on  numerical 
simulation  of  the  degradation  of  PEL  interaction  strength  due  to  these 
effects.  Results  of  the  gain  optimization  analysis  are  presented  in 
Section  III.  The  required  e-beam  properties  are  found  to  be  quite 
stringent  at  visible  wavelengths,  but  angular  canting  of  the  magnets  to 
provide  equal  two-plane  focusing  is  shown  to  enhance  emittance 
acceptance.  In  addition,  possible  tradeoffs  are  examined  which  relax  the 
e-beam  emittance  and  energy  spread  requirements,  but  produce  lower  gain 
per  unit  e-beam  current.  Implications  are  addressed  in  Section  IV. 


II.  Gain  Optimization 


A  general  approach  to  the  problem  of  optimization  of  the  PEL 
interaction  is  developed  in  this  section,  leading  to  specific  requirements 
for  the  wiggler  and  e-beam  properties.  The  gain-extraction  product  is  a 
useful  figure  of  merit  for  the  taper ed-wiggler  pel  interaction  strength 
since  gain  can  be  evenly  traded  for  extraction  (and  vice  versa)  by 
changing  the  wiggler  taper. 

A.  Gain-Extraction  Product 

The  gain-extraction  product,  { ,  has  been  shown^4^  to  be  given  by 
(cgs  units) 

i  -  (Osin?)2  >1x2  ,  [1] 


where 


1  h  w  In 

Xm  -  — 


(q  +  [1  +  q2]1/2) 


“  Xw  [l  +  a*]  (1  +  q2: 


6  is  the  fraction  of  electrons  trapped  in  the  ponderomotive  well,  sin?  is 
the  average  sine  of  the  phase  angle  for  trapped  electrons,  y  is  the 
electron  energy  in  units  of  the  rest  mass  me2,  i  is  the  electron  current, 
a  is  the  ratio  of  magnet  half-gap  h  to  the  1/e  photon  beam  amplitude 
radius  we  at  the  wiggler  entrance  and  exit,  aw  «  eB0Xw/23/2  tfmc?,  B0  is 
the  peak  B-field,  Xv  is  the  wiggler  wavelength,  q  -  1^/22^,  L*  is  the 
wiggler  length,  and  Zg  is  the  Rayleigh  range  of  the  photon  beam.  Equation 
[1]  does  not  include  degradation  of  the  gain-extraction  product  due  to 
energy  spread  and  emittance;  these  effects  will  be  introduced  later. 
Equation  [1]  is  derived  under  the  assumptions  of  small  fractional  change 
in  the  resonant  energy. 
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and  low  gain.  The  electron  dynamics  are  idealized  by  assuming  that 
trapped  electrons  behave  as  the  synchronous-phase  particle^5)  does,  so 
that,  in  essence  the  optimization  is  for  a  single  electron  traveling  along 
the  wiggler  axis. 

The  assumed  geometry  is  shown  in  Figure  1.  A  diffraction-limited 
photon  beam  is  focused  in  the  center  of  a  planar  wiggler.  A  Balbach^6^ 
magnet  configuration  is  used.  The  permanent  magnets  have  polarization 
vectors  oriented  as  indicated  by  the  arrows  in  Figure  1.  In  order  to 
properly  reflect  the  physical  limitations  of  developing  luge  magnetic 
fields  in  this  geometry,  we  write  av  in  terms  of  the  basic  magnet 
parameters  h,  Xw,  and  the  remanent  magnetization,  Br.  For  the  magnet  bar 
height  g  -  3Xw/8  and  no  gaps  between  neighboring  magnets,  the  calculated 
field  may  be  expressed  as^6) 

-2ffh/X 

av  -  1.07  x  io-«  Br  Xwe  w  .  [4] 

Vain g  a  typical  peak  value  of  Osin?  of  0.26  (based  on  numerical 
simulation  results),  a  conservative  clearance  factor  a  of  2,  and  Br  -  9000 
G  typical  of  SmCo5  magnets,  the  gain-extraction  product  is  now  a  function 
of  by,  Xw,  h,  X8,  ZR,  and  I.  Optimization  of  expression  Equation  [l] 
depends  on  which  parameters  are  fixed  and  which  are  varied  to  produce  the 
maximum.  The  case  of  X8,  L*,  and  I  fixed  is  considered  here.  These 
choices  allow  the  wiggler  length  to  be  held  to  a  practical  value  and 
assume  that  the  peak  current  available  is  limited .  It  is  to  be  noted  that 
when  optical  component  damage  is  an  issue,  ZR  may  be  a  more  useful 
independent  parameter  than  l^,  since  zR  partly  determines  the  beam  size  on 
the  resonator  mirrors.  Fixing  X8,  L^,  and  I  apparently  leaves  a  three- 
dimensional  surface  to  be  examined,  but  h  and  ZR  are  not  independent  of 
one  another  (due  to  the  fixed  clearance  factor  a)  and  the  resulting  two- 
dimensional  space  is  easily  analyzed  to  find  the  system  parameters 
yielding  the  maximum  gain-extraction  product. 
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The  e-bean  energy  spread  and  emittance  requirements  can  also  be 
determined  as  a  function  of  the  system  parameters.  The  allowable  energy 
spread  is  determined  by  the  requirement  that  it  cannot  exceed  the 
ponderomotive  bucket  full  height.  Similarly,  the  allowable  emittance  is 
set  by  the  combined  requirements  that  the  effective  energy  spread  due  to 
emittance  be  less  than  the  bucket  height  and  that  the  electron  and  photon 
beams  overlap  spatially. 

B.  Energy  Spread  Requirements 

For  a  tapered-wiggler,  the  energy  spread  must  be  less  than  the 
ponderomotive  bucket  height  full  width,  given  by^ 5  ^ 


f  2  X  e  a  Fry  I 
s  s  w  LrrJ 


Bucket 


it  l  +a 


where  es  **  eE0/21/2m0c2  is  the  normalized  rms  E-field,  F(yr)  *  cos fz  - 
(77/2  -  yr)sinyr,  and  parameter  fx  is  the  phase  angle  in  the  ponderomotive 
potential  well  for  the  electron  whose  energy  loss  rate  exactly  matches 
that  of  the  wiggler  resonant  energy.  A  phase  angle  of  about  40  degrees 
maximizes  the  product  of  trapping  fraction  and  bucket  deceleration  rate 
for  the  case  of  a  monoenergetic,  zero  emittance  beam.  At  photon  flux 
levels  too  low  to  result  in  trapping,  the  allowable  energy  spread 
corresponds  to  the  homogeneous  small-signal  linewidth.  Using  the 
linewidth  defined  by  Brau^ 7 )  for  the  linearly  tapered  wiggler  gives  an 
allowable  energy  spread  of 
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where  A>r  is  the  resonant  energy  change  of  the  wiggler  taper.  It  is  of 
interest  to  compare  the  energy  spread  requirements  at  small-signal  and 
saturated  flux  levels.  Since  saturation  (the  onset  of  trapping)  occurs 
for  an  E-field  value  of  roughly 
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the  bucket  full  width  at  saturation  is ,  for  ?r  m  40  ° , 


which  is  equal  to  the  equivalent  energy-width  due  to  the  small-signal 
linewidth.  Thus,  the  energy  spread  requirements  for  a  linearly  tapered 
wiggler  at  small-signal  and  saturated  conditions  are  identical.  At  flux 
levels  well  above  the  onset  of  saturation,  the  energy  spread  requirements 
relax  as  the  bucket  size  grows. 

C.  Emittance  Requirements 

The  allowable  emittance  in  the  focusing  plane  of  a  planar  wiggler  is 
often  set  by  the  requirement  that  the  effective  energy  spread  due  to 
emittance  be  less  than  the  bucket  height.  For  such  cases,  in  which 
spatial  overlap  is  not  the  limiting  factor,  determination  of  the  allowable 
emittance  is  based  on  the  principle  that  electrons  with  slightly  different 
trajectories  interact  in  the  same  way  as  electrons  with  identical 
trajectories  but  slightly  different  energies.  The  trajectory  difference 
is  of  course  related  to  the  emittance,  with  the  consequence  that  the 
allowable  emittance  can  be  directly  related  to  the  allowable  energy 
spread . 

The  effective  energy  spread  due  to  emittance  may  be  determined  by 
examining  the  variation  in  transverse  momentum  and  B-field  experienced  by 
an  electron  which  executes  the  betatron  orbit^ 8 )  of  maximum  amplitude . 

The  effective  energy  spread  full  width  is  found  to  be 
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where  k^  m  277/ Xw  ie  the  wavenumber  of  the  wiggler,  k^  is  the  wavenumber 
for  betatron  oscillations  resulting  from  the  distributed  focusing  along 
the  wiggler,  and  the  normalised  emittance  eN  is  defined  as  yvxQ  where  r 
is  the  radius  and  6  is  the  half-angle  at  any  beam  focus.  The  two  terms  in 
Equation  [9]  are  due  to  detuning  from  the  resonance  condition  by  the  B- 
field  gradient  and  by  the  trajectory  angle  associated  with  the  betatron 
orbit,  respectively.  The  minimum  energy  spread  occurs  for  a  focusing 
strength  kg  »  a^c^y  which  is  precisely  the  natural  value  for  a  planar 
wiggler.  In  that  case  the  detuning  from  the  resonant  condition  for  each 
electron  is  independent  of  axial  position.  Electron  trapping  will  be 
relatively  inefficient  when  the  effective  energy  spread  given  by  Equation 
[9]  exceeds  the  bucket  full  width.  For  other  focusing  strengths,  the 
energy  spread  given  in  Equation  [9]  is  actually  a  peak  value  achieved  only 
at  certain  points  in  the  betatron  orbit.  In  this  case  detrapping  will 
occur  if  the  energy  spread  exceeds  the  bucket  height  and  the  synchrotron 
wavelength  is  much  shorter  than  the  betatron  wavelength.  The  latter 
condition  is  marginally  satisfied  for  the  systems  of  interest.  The 
requirement  that  the  effective  energy  spread  be  less  than  the  well  depth 
will  hereafter  be  called  the  "bucket  constraint." 

In  some  cases,  the  allowable  emittance  ny  be  limited  by  the  need  to 
maintain  spatial  overlap  in  the  focusing  plane.  This  occurs  when  the 
allowable  e-beam  radius  based  on  the  bucket  constraint  exceeds  the  photon 
beam  waist  size.  The  e-beam  radius  given  by  the  distributed  focusing  in 
the  wiggler  is 


reb  ‘ 
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This  result  is  valid  when  the  electrons  are  optimally  focused  at  the 
wiggler  entrance,  in  which  cue  the  e-beam  radius  is  length  independent 
within  the  wiggler.  As  shown  later,  a  reasonable  condition  for  minimal 
gain  degradation  at  fixed  extraction  is  the  requirement  that  the  e-beam 


radius  be  less  than  the  1/e  photon  intensity  radius,  that  is,  r^  < 


w/21/2,  where  w  is  the  1/e  photon  amplitude  radius.  This  requirement 
leads  to  an  allowable  emittance  of 


rX  ZJc 


€N  < 


sR/S 


in: 


The  requirement  to  maintain  spatial  overlap  between  the  electron  and 
photon  beams  shall  be  called  the  "overlap  constraint.” 


In  the  free-expanding  plane,  spatial  overlap  is  generally  the 
limiting  constraint.  Proper  overlap  may  be  provided  by  matching  the 
photon  and  e-beam  envelope  shapes,  leading  to  the  requirement 


6N  < 


[12; 


Comparison  with  Equation  [11]  shows  that  the  emittance  requirement  for 
overlap  is  more  severe  in  the  free-expanding  plane  for  those  systems  which 
have  greater  than  unity.  This  constraint  may  be  relaxed  if  e-beam 

focusing  is  provided  in  the  nominally  free-expanding  plane. 


D.  Two-Plane  e-Beam  Focusing 


Two-plane  focusing  can  be  provided  in  a  planar  wiggler  either  by 
external  quadrupoles  or  by  angular  rotation ^ 9,10 )  of  the  wiggler  magnets. 
The  former  method  allows  a  readily  adjustable  focal  strength  but  requires 
precise  alignment  to  insure  that  the  wiggler  and  focusing  elements  are 
coaxial.  Otherwise  the  beam  will  be  steered  off-axis  and  a  betatron 
oscillation  excited.  The  latter  method  ensures  that  the  focusing 
properties  are  properly  aligned  and,  in  principle,  could  be  used  with 
hybrid  ( SmCo.  plus  steel )  wigglers  which  do  not  allow  linear  superposition 


«... I- 1 

•  -  • 


of  external  fields.  Two-plane  focusing  using  magnet  canting  has  recently 
been  experimentally  demonstrated  in  a  SmCo5  wiggler.^10^  When  focusing  is 
introduced  in  the  wiggle  plane,  the  focal  strength  in  the  nominally 
focusing  plane  is  reduced  according  to 


[13] 


where  k^x  and  k^y  are  the  wiggle  plane  and  nominally  focusing  plane 
betatron  wavenumbers,  respectively,  and  *  awkw/>  is  the  nominal 
betatron  wavenumber  in  the  absence  of  canting.  The  best  spatial  match 
with  a  cylindrically  symmetric  photon  beam  occurs  with  equal  two-plane 
focusing,  in  which  case 


a  k 
w  w 

2l/2y 


[14] 


E.  Confirmation  of  Emittance  and  Energy  Spread  Requirements 


Equations  [5],  [9],  [11],  and  [12]  define  energy  spread  and 
emittance  requirements.  The  actual  degradation  in  FEL  interaction 
strength  for  energy  spread  and  emittance  values  which  approaches  these 
limits  is  computed  in  this  section.  Electrons  which  are  detuned  to  near 
the  edge  of  the  bucket  or  experience  lower  E-fields  due  to  off-axis 
trajectories  experience  a  somewhat  weaker  interaction. 


These  effects  may  be  quantified  by  numerical  integration  of  the 
equations  for  electron  energy  loss  and  electron  phase  in  the  ponderomotive 
potential  of  the  FEL  interaction. 


d y_ 
dz 


"Vw  ,in  (f-** 


1  +  a 


V 

+  r* 


[15] 


The  effects  of  emittance  are  included  in  a  two-step  process.  First  the 
transverse  angles  0X  and  By  of  the  electron  orbits,  the  optical  E-fields, 
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and  the  wiggler  B-fields  experienced  by  each  electron  are  computed  as  a 
function  of  axial  position,  including  the  off-axis  notion  but  ignoring  the 
snail  effect  of  the  PEL  interaction  on  the  trajectory.  Then  the  energy 
loss  is  determined  in  a  one-dimensional  integration  of  the  coupled 
equations  Equation  [15]  using  the  previously  computed  angles  and  E  and  B- 
fields  for  each  electron.  The  emittance  phase  space  is  assumed  to  be 
uniformly  filled.  The  electron  and  optical  beams  are  assumed  to  be 
optimally  focused  and  coaligned  and  the  optical  beam  is  assumed  to  be 
diffraction-limited,  in  which  case,  in  the  low-gain  approximation,  the  E- 
field  amplitude  and  phase  varies  spatially  according  to*11) 


,  v  o  [  wo  1  -(r/w[z])2 

Vr's>  ■  ••  IsiirJ* 

*(T.Z)  -  -  [jiy]  (VZ*)  . 


tie; 


where  e°  is  the  normalized  on-axin  E-field  at  the  waist,  w(z)  •  wQ  (l  + 
[z/ZR]i  )1/2,  wc  is  the  1/e  photon  amplitude  radius  at  the  waist,  and  z  is 
measured  from  the  waist  location.  The  spatial  B-field  variation  is 
included  by  assuming  constant  ^  with  aw  tapered  for  constant  resonant 
phase 

o 

V  y,z)  -  av  C  cosh  V  .  U7; 


where 


C(*)  -  1  - 


Aa 


In 

-q  +  (i  +  q2)1/2 _ 1 

,/z^  ♦  (l  *  [Z/Z*)  a) 

-T/Ij 

In 

,  1/2' 

-q  +  (i  +  q2 ) 

q  +  (1  +  q2 )  ' 

[is; 


a£  is  the  aw  value  at  the  entrance,  and  Aaw  is  the  change  in  a^,  along  the 
wiggler  length. 


Numerical  calculations  have  been  completed  for  the  preliminary  point 
design  of  a  0.5  /m  oscillator  experiment  given  in  Table  1.  This  design 
incorporates  a  SmCo5  wiggler  of  5  m  length  and  an  e-beam  energy  of  120 
MeV.  The  wiggler  design  assumes  two-plane  e-beam  focusing  provided  by  the 
wiggler  magnets . 

Results  of  the  numerical  electron  tracking  calculations  are  given  in 
Figure  2.  The  effect  of  emittance  is  to  require  an  overall  higher  photon 
power  to  achieve  a  given  level  of  extraction.  Since  overlap  is  the  more 
severe  emittance  constraint  (see  Table  1),  the  higher  photon  flux 
compensates  for  the  relatively  low  E-fields  experienced  by  most  electrons 
due  to  their  off-axis  trajectories.  While  the  additional  E-field  required 
to  achieve  extraction  at  higher  emittance  is  modest,  the  associated 
increase  in  photon  power  may  be  significant.  About  50  percent  gain 
degradation  is  found  at  the  overlap  constraint  of  r^  »  wQ/21/2  defined 
previously.  A  similar  gain  degradation  factor  is  obtained  for  an  energy 
spread  which  just  matches  the  full  bucket  height.  These  results  show  that 
the  previously  defined  energy  spread  and  emittance  requirements  correspond 
to  defining  the  acceptance  as  the  value  which  degrades  the  FEL  interaction 
strength  by  one-half. 


III.  Results  and  Discussion 

The  optimization  analysis  developed  in  the  previous  section  is  used 
to  define  the  parameters  for  systems  optimized  for  peak  gain  at  fixed 
extraction.  Energy  spread  and  emittance  acceptance  values  are  identified 
for  these  systems.  Later,  tradeoffs  which  lower  the  gain  but  enhance 
acceptance  are  examined. 

A.  Gain-Optimized  Systems 

The  calculated  optimum  gain  and  corresponding  e-beam  energy  are  given 
in  Figure  3.  The  gain  values  are  in  units  of  percent  optical  power 
increase  per  asipere  of  e-beam  current.  They  apply  to  5  percent  extraction, 
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Table  l 


I 

PRELIMINARY  POINT  DESIGN  FOR  0.5 
OSCILLATOR  EXPERIMENT 


E-Beam  Energy 

120  MeV 

Peak  Current 

100  A 

Viggler  Length 

5  m 

Taper  (Ayr/>r) 

12% 

Wiggler  Wavelength 

2.18  cm 

Peak  Magnetic  Field 

8.7  kG 

a 

w 

1.25 

Full  Gap 

0. 36  cm 

Rayleigh  Range 

2.2  m 

Instantaneous  Photon  Power  at  5%  Extraction 

3.3  GW 

Normalized  Emittance  Acceptance 

Based  on  Overlap 

0.014  cm-rad 

Normalized  Emittance  Acceptance 

Based  on  Effective  Energy  Spread3 

0.048  cm-rad 

Energy  Spread  Acceptance3 

1.5% 

.  b 

Single-Pass  Gain  at  5%  Extraction 

18% 

a 

Assumes  equal  two-plane  focusing  provided  by  wiggler 
and  50  percent  gain  loss  at  fixed  extraction. 

» 

Assumes  €w  ■  0.01  cm-rad,  Energy  Spread  «  1  percent. 
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HD-M.44  952 
UNCLASSIFIED 


TAPERED-H I QOLER  FREE-ELECTRON  LASER  OSCILLATOR  PROORAH  2/2 
(U>  NATHENATICAL  SCIENCES  NORTHWEST  INC  BELLEVUE  HA 
J  SLATER  ET  AL.  HAY  84  HSNH-C.  11.  19S.  81 

AF0SR-TR-84-8E23  F49628-81-C-8R79  F/Q  29/5  NL 


at  5%  Extraction 


Figure  3.  Maxima  Single-Pax*  Gain  and  optimum  E-beam  Energy. 
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-Beam  Energy  (MeV) 


but  can  be  trivially  scaled  to  other  extraction  values  since  the  gain- 
extraction  product  is  constant.  Comparison  with  the  numerical  tracking 
code  described  previously  shows  that  these  gains  are  generally  accurate  to 
within  30  percent.  The  curves  clearly  indicate  that  higher  gains  may  be 
obtained  with  longer  wigglers.  This  results  primarily  £rom  the  term 

+  «£>  in  Equation  [2],  which  increases  in  value  as  the  bore  size 
h  increases.  The  gain  is  not  a  strong  function  of  the  photon  wavelength. 
This  is  due  to  our  assumption  of  constant  current,  which  means  that  the 
weaker  interaction  at  shorter  wavelengths  is  partially  offset  by  the  larger 
e-beam  power  at  higher  y  values.  The  associated  optical  powers,  P8,  for  5 
percent  extraction  are  given  in  Figure  4.  These  powers  can  be  scaled  to 
extraction  values,  7 j,  other  than  5  percent  by  noting  that  tj3/P8  is 
constant . 

optimum  values  of  the  dependent  parameters  for  the  1  micron  photon 
wavelength  case  are  shown  in  Figure  5.  Parameters  L^/Zp  and  aw  remain 
roughly  constant  while  the  ratio  h/Xw  varies  significantly  with  wiggler 
length.  It  is  interesting  to  note  that  the  optimum  value  is  less  than 
unity  for  the  conditions  examined  and  for  this  particular  choice  of  the 
independent  parameters .  An  optimum  aw  value  of  1.0  is  frequently  reported 
in  the  FEL  literature,  and  the  latter  value  doeB  result  from  this 
optimization  when  Xw  is  taken  to  be  independent  (i.e.,  fixed).  Certainly 
one  can  pick  Xw  to  be  a  system  constraint  if  so  desired,  but  equivalent 
wiggler  performance  will  then  require  slightly  longer  wigglers  than  for  the 
case  where  is  constrained  and  Xw  is  optimized. 

Calculated  values  of  the  energy  spread  requirements  for  gain- 
optimized  systems  are  shown  in  Figure  6.  These  curves  apply  to  the  5 
percent  extraction  case,  but  may  be  scaled  to  other  values  by  noting  that 
the  bucket  height  varies  proportional  to  T?1/3.  The  bucket  depth 
dependence  on  photon  wavelength  is  weak  because  the  XBe8  product  in 
Equation  [5]  does  not  vary  substantially.  These  bucket  widths  have  been 
calculated  for  the  field  at  the  entrance  or  exit  of  the  wiggler.  From  the 
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Higgler  Length  (m) 
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Figure  6 


Energy  Spread  Requirement  at  Optimum  Gain  Conditions 
and  5  Percent  Extraction  Based  on  Equating  Maximum 
Energy  Spread  to  Bucket  Height. 


Ly/Zji  values  of  Figure  5,  the  bucket  at  the  wiggler  center  is  typically  50 
percent  larger. 

The  emittance  values  which  just  fill  the  bucket  and  Which  just  meet 
the  overlap  requirement  are  shown  in  Figure  7.  Clearly  the  emittance 
requirement  for  overlap  is  the  more  severe  constraint  and  becomes 
particularly  serious  at  shorter  photon  wavelengths .  This  constraint  may  be 
relaxed  somewhat  if  e-beam  focusing  is  provided  in  the  nominally  free- 
expanding  plane. 

B.  Two-Plane  e-Beam  Focusing 

The  allowable  emittance  for  gain  optimised  systems  with  equal  two- 
plane  focusing  is  shown  in  Figure  8.  As  with  single-plane  focusing,  the 
allowed  emittance  is  generally  limited  by  spatial  overlap  constraints.  For 
short  wiggler  lengths,  the  overlap  obtained  with  two-plane  focusing  and  an 
axially  independent  beam  size  is  worse  than  the  overlap  Which  can  be 
achieved  with  no  wiggler  focus  and  a  slightly  converging  input  beam.  In 
this  case  the  emittance  requirement  shown  is  identical  to  that  in  Figure  7 
for  single-plane  focusing.  For  Xg  -  10  (u n,  the  allowed  emittance  is 
limited  by  the  bucket  constraint  for  wiggler  lengths  greater  than  8 
meters.  For  the  longer  wiggler  lengths,  the  emittance  acceptance  is 
considerably  enhanced  relative  to  the  single-plane  focusing  case.  Since 
the  actual  improvement  in  useful  current  scales  as  c3,  two-plane  focusing 
represents  a  significant  advantage. 

The  overlap  problem  may  be  further  affected  by  additional  external 
focusing  along  the  wiggler  length.  Such  additional  focusing  may  be  useful 
whenever  the  tural  focus  of  the  wiggler  produces  an  effective  energy 
spread  that  does  not  exceed  the  bucket  depth.  In  Figure  9,  the  emittance 
requirements  based  on  the  bucket  (Equation  [9])  and  overlap  (Equation  [11]) 
constraints  are  shown  as  a  function  of  focusing  strength  by  the  clashed  and 
solid  lines,  respectively.  This  example  case  has  XB  -  0.5  im,  5  percent 
extraction,  and  1^  ■  5  m.  Since  the  emittance  acceptance  of  the  wiggler  is 
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Figure  7.  Normalized  Emittance  to  Fill  Bucket  in  Focusing  Plane 
and  for  Spatial  Overlap  in  Free-Expanding  Plane. 
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Higgler  Length  (m) 


Figure  8 


Emittance  Acceptance  of  Gain-0ptimi2ed  Systems  with 
Equal  Two-Plane  Focusing  Provided  by  Higgler.  Bucket 
constraint  based  on  5  percent  extraction. 


Normalized  Emittance  Acceptance  lyurQ,  cm-rad) 


the  minimum  of  the  two  constraints,  the  optimnn  operating  point  is  expected 
to  be  at  their  intersection.  Strong  external  focusing  appears  to  provide  a 
promising  mans  of  emittance  acceptance  enhancement. 


C.  Optimization  for  High  nti.ttanoe 

In  the  previous  section  systems  were  optimized  for  highest  gain  under 
the  assumption  of  zero  emittance,  and  the  emittance  acceptance  was  then 
computed .  For  nonzero  emittance  values,  the  wiggler  parameters  specified 
do  not  provide  the  highest  possible  gain.  An  optimization  for  larger 
emittance  values  is  now  considered  for  systems  with  equal  two-plane 
focusing  provided  by  the  wiggler.  For  simplicity  the  analysis  is  tailo- 
to  cases  where  the  equivalent  energy  spread  due  to  emittance  is 
unimportant,  as  is  the  case  for  the  following  0.5  im  calculations. 


The  function  to  be  optimized  is  equivalent  to  that  considered 
previously  (Equation  [1]),  except  that  new  factors  are  added  to  reflect  the 
degradation  of  gain  at  fixed  extraction  with  increasing  emittance  and 
energy  spread.  For  a  given  fractional  energy  spread  AE  -  Ay/y  and 
emittance  eN,  these  factors  depend  on  the  other  system  parameters  and  are 
given  by 


[19] 


[20] 


where  h  is  the  bucket  height  defined  by  Equation  [5]  and  e0  is  the 
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normalized  emittamce  acceptance  based  on  the  overly  constraint  given  in 

Equation  [11] .  The  functional  dependences  for  lower  energy  spread  and 

omittance  values  ( Ae/H  <  (2/3 )1/a ,  eN/eo  <  4/3)  axe  fits  to  the  numerical 

l  /2 

results  shown  in  Figure  2 .  For  large  energy  spread  ( AE/H  >  ( 2/3 )  '  ),  the 
fraction  of  electrons  in  the  bucket  is  proportional  to  1/AE,  and  the 
proportionality  constant  is  chosen  to  connect  smoothly  to  the  curve  in 
Figure  2(b).  This  scaling  properly  refers  to  cases  where  energy  spread  in 


excess  of  (2/3) 


V2 


H  is  filtered  out  upstream  of  the  wiggler.  For  large 


emittance  ( CN/€Q  >  4/3)  the  fraction  of  electrons  within  the  optical  beam 
is  proportional  to  eN-2,  and  the  proportionality  constant  is  chosen  to 
connect  smoothly  to  the  curve  in  Figure  2(a). 


The  new  function  to  be  optimized  is  then 

,  12877 2 e  , ,  .  2  2  . 

«  “  (OBint)  FIX  f6Efe 


[21] 


As  before,  this  is  an  expression  for  the  gain  at  fixed  extraction.  The  two 
additional  parameters  AE  and  cN  are  taken  to  be  fixed,  while,  as  before, 

Xw,  h,  and  Zp  are  varied  to  find  the  optimum.  The  results  of  this 
analysis  for  Xfl  -  0.5  im  are  shown  in  Figure  10.  For  wiggler  lengths  over 
4  meters,  the  system  parameters  adjust  to  accept  normalized  emittance 
values  of  up  to  0.01  cm-rad  with  no  more  than  50  percent  gain  loss.  Energy 
spreads  of  up  to  1  percent  result  in  no  more  than  50  percent  gain  loss  for 
wiggler  lengths  of  under  7  meters.  Intermediate  wiggler  lengths  of  4  to  7 
meters  are  leant  sensitive  to  the  combined  effects  of  emittance  and  energy 
spread . 


Figure  li  shows  how  system  parameters  change  to  accept  larger 
emittance.  To  provide  proper  overlap  as  the  emittance  is  increased,  the 
photon  beam  size  must  be  increased  beyond  the  optimal  zero-emittance  size. 
In  addition,  the  e-beam  energy  is  increased  because  the  equivalent  energy 
spread  decreases  with  y.  The  photon  beam  size  increase  causes  gain  loss 
in  two  ways.  First  the  total  photon  power  increases  proportionally  to  the 
beam  area  so  that  less  energy  is  extracted  per  unit  photon  energy 
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Figure  10.  Maximum  gain  at  As  *  0.5  ym  and  5  Percent  Extraction 
for  Various  Assumed  Levels  of  Bnittance  and  Energy 
Spread. 
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Fraction  of  Value  at  Optimum  Gain  Condition 


Relative  Quittance  Acceptance 


Figure  11.  Variation  of  System  Parameters  to  Accept  Larger  Emittance 
While  Maintaining  as  High  a  Gain  as  Possible.  L  »  5  m, 

*  *  0.5  pm,  5  percent  extraction,  AE  =  0.  w 


invested.  Second,  the  wiggler  bore  size  must  also  increase,  thereby 
reducing  the  B-field  amplitude  and  the  interaction  strength,  we  note  that 
the  magnet  gap  actually  decreases  at  first  as  the  photon  radius  increases. 
This  is  because  the  gap  is  a  fixed  multiple  of  the  photon  beam  size  at  the 
wiggler  ends,  and  as  moves  from  the  optimum  value  of  about  4  to 

smaller  values,  the  end  size  decreases  to  a  minimum  value  When  the  wiggler 
is  confocal  (1^  -  22^).  Not  surprisingly,  the  curves  show  roughly  that  the 
emittance  acceptance  can  be  improved  fourfold  by  doubling  the  photon  beam 
radius  at  the  expense  of  a  factor  of  2  gain  loss.  Ihis  benefit  to  cost 
ratio  becomes  less  favorable  for  further  increases  in  the  photon  beam 
radius  and  magnet  gap. 

It  is  of  interest  to  compare  the  falloff  in  gain  with  emittance  for 
various  designs.  The  approximate  falloff  with  increasing  emittance  is 
given  directly  by  Equation  [20],  Figure  12  compares  a  set  of  these  rolloff 
curves  for  various  5  meter  wiggler  designs  operating  at  0.5  um  and  5 
percent  extraction.  These  wigglers  differ  in  gap  and  magnetic  field 
wavelength  so  as  to  provide  different  emittance  acceptance  values.  It  is 
evident  in  the  figure  that  over  the  range  shown  the  gain  at  zero  emittance 
is  falling  in  rough  proportion  to  the  square  root  of  the  increase  in 
emittance  acceptance,  assuming  the  acceptance  is  defined  as  the  point  at 
which  the  gain  falls  by  one-half. 


IV.  Implications 

New  subharmonic  bunching  linac  injectors  recently  developed  at  the 
Boeing  Aerospace  Company  (BAc/12^  and  at  the  Stanford  Linear  Collider 
(SLC)^13)  have  each  demonstrated  simultaneous  achievement  of  high  peak 
current  (>100A)  and  low  normalized  emittance  (<0.02  cm-rad).  These 
achievements  can  be  compared  to  the  well-known  Lawson-Penner 
relation^ 1-3 ' 9  ^ 


2 

<I>  [kA]  -  e  [cm-rad]  , 

N 


[22] 
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Figure  12.  Gain  Versus  Quittance  for  Systems  Optimized  for 

Various  Levels  of  Quittance.  L  =  5  m,  \  «  0.5  pm 

5  percent  extraction,  AE  *  0.  w  8 


which  is  observed  to  approximately  represent  the  relationship  between  the 
time-averaged  current  <X>  and  the  normalised  emittance  of  a  wide  range  of 
accelerators  working  at  average  currents  from  milliamperes  to  kiloamperee. 
RF  accelerators  generally  produce  some  degree  of  bunching  while  the 
electrons  are  of  relatively  low  energy,  so  the  peak  and  average  current 
values  differ,  lhe  microscopic  duty  cycle  of  the  high  peak  current 
injectors  mentioned  previously  is  made  especially  small  (<0.01)  through 
the  use  of  subharmonic  bunching,  yielding  high  peak  currents  with 
comparatively  low  emittance. 

The  implications  of  this  accelerator  technology  development  on  PEL 
performance  is  summarized  in  Figure  13.  Here  the  single-pass  gain  at  5 
percent  extraction  is  plotted  as  a  function  of  peak  current  and  normalised 
emittance  for  the  preliminary  point  design  of  a  visible  oscillator  given  in 
Table  1.  that  design  was  optimised  for  the  8AC  accelerator,  and  as  shown 
in  the  figure,  about  20  percent  single-pass  gain  is  possible  at  5  percent 
extraction,  included  for  comparison  are  the  SVC  accelerator  and  the 
Stanford  superconducting  linac/14)  which  w as  used  for  the  first  FEL 
experiment/15^  All  three  accelerators  are  characterised  by  roughly  the 
same  current  densities  in  emittance  phase  space,  (lp/ena),  although  the 
peak  currents  differ  by  two  orders  of  magnitude.  State-of-the-art  high 
peak  current  RF  linacs  appear  to  satisfy  the  e-beam  requirements  for  a  high 
extraction  oscillator  at  visible  wavelengths. 


Normalized  Emittance  (yirrS,  cm-rad) 


63  07673 


Figure  13.  Recent  Improvements  in  Emittance  and  Peak  Current  Achieved 
by  Linear  Accelerators  Improve  the  Prospects  for  Visible 
Oscillation  with  High  Efficiency.  Gain  values  at  0.5  vim 
assure  L  “5m  and  5  percent  extraction. 
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APPENDIX  C 


Mode  Structure  of  a  Tapered-wiggler 
Free-Elect ron  Laser  Stable  Oscillator 


The  tapered-wiggler  free-electron  laser  (PEL)  is  of  current  interest 
for  its  potential  as  an  efficient,  tunable  laser.  The  basic  idea  is  a 
variation  of  the  PEL  concept  demonstrated  in  the  experiment  of  Deacon,  et 
al.[l],  in  which  it  was  shown  that  a  small  fraction  of  electron-beam 
kinetic  energy  could  be  converted  to  optical  energy  in  a  single  pass 
through  a  wiggler  magnet.  The  tapered-wiggler  concept[2]  involves 
variation,  or  tapering,  of  the  wiggler  parameters  as  a  function  of  axial 
position  in  such  a  manner  as  to  maintain  a  resonant  interaction  as  the 
electrons  decelerate.  This  tapering  allows  increased  kinetic  energy 
extraction  per  pass,  and  may  lead  to  high  efficiency  systems  through 
reduction  of  e-beam  recirculation  or  energy  recovery  requirements. 

Several  verification  experiments  of  the  basic  concept  are  in 
progress[3].  Initial  experiments  are  expected  to  achieve  on.'y  low  photon 
gain,  and  in  a  sense  provide  validation  only  of  the  relatively  simple 
theory  describing  electron  deceleration  in  predetermined  photon  and 
wiggler  fields.  The  process  of  self-cons istently  determining  electron  and 
photon  behavior  in  high  gain  amplifiers  or  oscillators  is  significantly 
more  complex,  and  has  recently  been  the  subject  of  considerable  analysis . 
With  respect  to  the  photon  beam  parameters,  such  analysis  has  generally 
fallen  along  one  of  two  lines,  study  of  either  the  axial  profile  of  the 
beam  When  short  pulses  are  of  interest[4],  and  study  of  the  transverse 
structure[5,6] .  The  theoretical  analysis  presented  here  pertains  to  the 
transverse  optical  field  structure  of  tapered-wiggler  oscillators. 

Interest  in  the  oscillator  transverse  field  structure  stems  largely 
from  the  nonlinear  interaction  of  the  gain  media  and  optical  field. 
Geometries  optimized  for  maximum  electron-photon  interaction  will 
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generally  involve  a  long,  mall  diameter  e-beam  centered  in  a  larger 
photon  beam.  The  resulting  amplitude  and  phase  modification  of  the  photon 
beam,  being  radially  nonuniform,  leads  to  generation  of  optical  cavity 
SK>des  which  should  be  accounted  for  in  system  design,  and  also  leads  to  an 
output  beam  of  lower  quality  than  the  ideal  diffraction  limit.  It  is  the 
quantification  of  these  effects  that  is  pursued  here. 

The  method  of  calculation  is  described  in  Sec.  II .  Phase  and 
amplitude  profiles  of  the  photon  field  are  determined  in  a  plane 
perpendicular  to  the  cavity  axis,  this  plane  being  numerically  propagated 
back  and  forth  along  the  axis  using  a  paraxial  wave  equation,  section  III 
describes  the  application  of  this  formalism  to  the  PEL  interaction  and 
relevant  geometries.  The  complexity  of  the  PEL  interaction  is  reduced  by 
means  of  a  resonant-phase  approximation,  thus  providing  a  simple 
relationship  between  the  phase  of  the  stimulating  and  stimulated  electric 
fields.  While  the  calculations  presented  have  relatively  simple  modal 
structure,  the  analysis  is  directly  applicable  to  more  complicated 
structure.  The  resonant-phase  approximation  does  limit  the  applicability 
of  the  model  to  fully  saturated  photon  intensity  levels.  The  optical 
cavity  is  taken  to  be  injected  from  an  external  source  at  the  full 
saturated  intensity,  and  the  subsequent  development  of  mode  structure  is 
the  subject  of  primary  interest.  Actual  tapered-wiggler  PEL  oscillators 
may  be  self  started  from  spontaneous  emission  or  require  injection  at  or 
below  saturation  intensity.  Our  analysis  method  precludes  modeling 
transverse  mode  structure  effects  in  the  small-signal  start-up  regime,  but 
it  provides  a  useful  way  to  determine  the  steady  mode  structure  in  the 
saturated,  trapped-particle  regime  of  tapered-wiggler  PEL  oscillators. 
While  the  details  of  the  mode  structure  during  the  approach  to  steady 
state  will  be  dependent  on  the  start-up  method,  the  steady-state  structure 
should  be  insensitive  to  the  start-up  method.  The  properties  of  the 
transverse  field  are  analyzed  by  means  of  decomposition  into 
Gaussian-Laguerre  modes,  providing  a  convenient  measure  of  beam  properties 
in  terms  of  partition  of  the  total  power  into  modes.  The  mode  structure 


depends  strongly  on  the  degree  of  aperturing  in  the  optical  cavity,  this 
resulting  frosi  finite  sise  optical  elements.  Aperturing  toy  the  viggler 
aagnets  will  toe  particularly  important  in  geometries  optimised  for  FEL 
interaction  strength [7 ] .  Presence  of  the  wiggler  has  been  accounted  for 
toy  truncation  of  the  beam  at  positions  of  the  beam  corresponding  to  the 
wiggler  entrance  and  exit  locations.  Results  of  the  atudy  are  presented 
in  Sec.  XV  and  sumarized  in  Sec.  v. 

XX.  CALCUIATIOHAL  METHODS 

The  free-electron  laser  oscillator  is  modeled  using  an  axisymmetric 
physical  optics  code,  within  the  wiggler,  the  wave  front  is  propagated 
numerically  using  a  finite-difference  solution  of  the  paraxial  wave 
equation  [8].  For  propagation  between  the  wiggler  and  distant  mirrors, 
the  Huygens-Fresnel  integral  [ 9]  is  used.  The  computer  model  is  used  to 
find  the  transverse  mode  structure  and  far-field  characteristics, 
accounting  for  diffraction,  nonuniform  gain  media,  refraction,  and 
arbitrary  mirror  configurations  and  reflectivities. 

A  linearly  polarized  wave  amplitude  may  be  expressed  as  a  scalar 
function  of  position  and  time 

u(x,y,z,t)  -  U(x,y,z)cos(wt  +  «(x,y,z)),  [1] 

where  U,  u,  and  $  are  the  amplitude,  angular  frequency,  and  phase  of  the 
wave,  respectively.  Using  phasor  notation,  the  explicit  time  dependence 
may  be  removed 

U  -  U(x,y,z )exp(-i$(x,y,z) ),  [2] 

where  we  use  the  superscript  to  denote  a  complex  number,  or  phasor. 

The  original  time-dependent  amplitude  may  toe  recovered  toy  the 
transformation 

u(x,y,z,t )  -  Re[U(x,y,z)exp( -iwt )] . 


[3) 


Trie  complex  wave  amplitude  propagates  through  free  space  according  to 
the  time- independent  wave  equation 


2  ~  2~ 

V  u  +  X  U  -  0, 


[4] 


where  k  •  2 n/\  is  the  laser  wave  number.  For  light  traveling  nearly 
parallel  to  the  z-axis,  one  can  remove  most  of  the  rapid  phase  variation 

A 

by  defining  a  new  function  y  by 


U(x,y,z)  -  y(x,y,z)exp(-i)cz). 


[5] 


The  function  y  represents  the  difference  between  the  actual  wave  front  and 
a  uniform  plane  wave,  and  thus  is  a  much  slower  varying  function  of 
position,  mis  function  satisfies 


2  -  ay 

V  y  -  2Dc  — ■  -  o. 
dz 


[6] 


the  so-called  "paraxial  wave  equation".  We  solve  this  equation 


numerically  in  axisymnetric  cylindrical  geometry  under  the  assumption  that 


a  2  2 

y  varies  so  slowly  with  z  that  its  second  derivative,  d  y/flz  ,  may  be 


neglected, 

1  a  f  ay]  ay  „ 

r  dr  lr  drl  dz 


[7] 


This  free-space  propagation  equation  is  solved  implicitly  using  a 
centered-dif ference  technique  involving  inversion  of  a  tridiagonal 
matrix[10].  This  method  is  stable  independent  of  step  size,  thus  allowing 
the  step  size  Az  to  be  chosen  on  the  basis  of  desired  accuracy,  rather 
than  on  stability  considerations.  The  technique  is  second-order  accurate 
and  is  found  to  achieve  excellent  energy  conservation. 


A  disadvantage  of  the  finite-difference  propagation  technique  is  the 
inability  to  accurately  treat  sharp  aperture  edges.  At  an  aperture,  the 
complex  amplitude  is  discontinuous,  so  that  high  order  radial  derivatives 
are  undefined.  Truncation  errors  due  to  these  terms  cannot  be  made 
arbitrarily  small  by  decreasing  Ar.  This  difficulty  is  overcome  by  using 


smoothly  tapered  aperture  edges[ 8 ] ,  allowing  accurate  numerical 
calculation  of  the  diffraction  pattern.  The  distinction  between  smooth 
and  hard  apertures  is  unimportant  to  this  application  because  differences 
in  the  fields  are  found  only  at  locations  very  close  to  the  aperture[ll] . 
The  apertures  are  smoothed  by  using  a  Gaussian  profile  at  the  edge.  Por 
example,  the  mirror  reflectivity  is 


r  <  r 


[8] 


R(r) 


Roexp[(r-to)I/Tll  r>V 


where  rQ  is  the  radius  of  the  portion  of  the  mirror  with  uniform 
reflectivity,  Rq,  and  r  is  the  truncation  distance.  We  choose  rc  and  r  so 
that  the  actual  diameter,  23^,  of  the  mirror  being  modeled  lies  at  the 
half  power  points  of  the  Gaussian  edges. 


In  regions  of  the  problem  where  the  PEL  gain  medium  is  active,  the 
propagation  equation  contains  an  additional  term  describing  the  PEL 
interaction 


11 L  fr  if]  _2ik  [air  _ 

r  3r  lr  dr  Uz  AtC  r 


) 


-i  (7T/2  -0)1 


1*1 


[9] 


Here  parameters  A  and  0(r,z)  describe  the  E-field  generated  by  the  gain 
medium.  The  parameter  A  is  the  amplitude  of  the  E-field  generated 
on-axis,  per  unit  distance  along  the  wiggler  (i.e.,  am  interaction  of 
length  dz  generates  a  field  of  amplitude  A  dz).  A  is  proportional  to  the 
density  of  trapped  electrons.  Since  the  number  of  electrons  in  the  bucket 
is  fixed,  A  is  independent  of  axial  position.  The  variation  in  gain  along 
the  wiggler  length  is  due  to  the  variation  of  the  phase  0  of  the 
stimulated  field  (the  field  of  amplitude  A  dz)  with  respect  to  the 

/v 

stimulating  field,  *.  Phases  of  490  and  -90  degrees  correspond  to  gain 
and  absorption,  respectively,  with  no  phase  shift.  Phases  of  0  and  180 
degrees  correspond  to  phase-front  lag  and  advance,  respectively,  with  no 
gain.  The  dimensionless  f(r)  is  a  spatial  fora  factor  describing  the 
radial  electron  density  profile. 


Til*  gain  ten  la  included  explicitly  in  the  wave-front  propagation 
algorithm  using  an  iterative  procedure  to  achieve  centered  differencing 
and  second-order  accuracy.  The  propagation  algorithm  is  solved  on  a 
uniform  two-dimensional  grid  with  about  200  steps  in  the  radial  direction 
and  100  steps  in  the  axial  direction.  The  problem  is  solved  by  specifying 
an  initial  wave  front  at  a  particular  axial  position,  and  then  using  the 
propagation  algorithm  to  advance  the  wave  front  forward  in  s 
simultaneously  at  each  radial  position.  When  necessary,  due  to  a  more 
rapidly  varying  complex  electric  field,  the  axial  grid  is  subdivided 
locally. 

Several  checks  confirm  the  accuracy  of  the  finite-difference 
propagation  algorithm.  Control  calculations  have  been  made  using 
different  grid  sizes.  These  calculations  verify  that  the  results  are 
independent  of  grid  site,  for  the  grid  resolution  used  in  this  study. 
Furthermore,  the  beam  power  is  conserved  to  within  0.05  percent  when  a 
combination  of  TEMqq  and  TEM10  modes  is  propagated  for  35  round  trips 
around  an  unape  rtured  confocal  resonator. 

SIX .  APPLICATION  TO  PEL  GEOMETRIES 

The  PEL  cavity  geometry  is  shown  in  Pig.  1. 


Figure  1.  Schematic  of  axisymmetric  FEL  oscillator  geometry.  Racial 
dimension  is  greatly  exaggerated  for  clarity.  Zr  ■  Rw*/X. 


A  plane-polarized  optical  bean  is  circulated  in  a  symmetric  stable 
resonator.  The  wiggler  length  is  equal  to  two  Rayleigh  ranges  of  the 
optical  bean.  This  geonetry  produces  a  shallow  waist  at  the  wiggler 
center,  and  nearly  optimizes  the  electron-photon  energy  exchange.  The 
majority  of  calculations  presented  in  this  paper  pertain  to  confocal 
cavities,  in  which  case  the  mirrors  are  located  precisely  at  the  ends  of 
the  wiggler.  In  practice,  the  mirrors  may  be  separated  by  more  than  the 
wiggler  length  to  provide  e-bean  access,  but  the  basic  aspects  of  the 
problem  remain  the  same.  The  cavity  mirrors  have  a  radius  of  a^.  The 
presence  of  the  wiggler  magnets  is  accounted  for  by  including  additional 
apertures  of  radius  a^  at  each  end  of  the  wiggler.  For  the  following 
analysis,  the  mirror  on  the  right  is  taXen  to  be  a  partially  reflective 
output  coupler,  with  radially  independent  reflectivity.  The  output 
coupling  is  chosen  to  be  approximately  equal  to  the  FEL  gain  at 
saturation.  Low  to  moderate  gain,  less  than  50  percent  per  pass,  is 
considered  in  all  cases.  The  electron  beam  is  assumed  to  have  a  Gaussian 
radial  electron  density  profile 

f ( r )  ■  exp[-r2/r^J,  [10] 

independent  of  axial  position.  In  all  calculations  presented,  the  e-beam 
is  assumed  to  be  fairly  small  in  diameter  compared  to  the  photon  beam. 

Calculation  of  the  interaction  with  the  gain  medium  is  based  on  the 
synchronous-particle  approximation.  That  is,  all  interacting  electrons 
are  assumed  to  be  trapped  in  the  ponderomot ive  potential  well,  and  it  is 
further  assumed  that  all  such  electrons  are  located  in  the  well  at  the 
stable-phase  point[2].  This  in  turn  determines  the  relative  phase  between 
the  locally  generated  field  and  locally  applied  field,  with  the  result 
that  6  in  Eq.  [9]  can  now  be  replaced  with  the  stable-phase  angle.  This 
approximation  is  useful  since  trapped  electrons  oscillate  about  the 
synchronous-phase  point,  but  its  use  precludes  modeling  effects  associated 
with  the  pendulum  motion  of  trapped  particles. 


_  • 
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The  stable-phase  angle,  ®,  is  a  function  of  the  local  electric  field 
and  local  viggler  taper  d?r/dz,  where  Yr  is  the  resonant  energy.  The 
wiggler  taper  must  be  specified  as  an  input  parameter.  For  this  analysis, 
the  taper  is  chosen  so  that,  in  the  presence  of  a  TEMq0  photon  bean,  the 
phase  would  be  independent  of  axial  position.  In  the  presence  of  the  FEL 
interaction,  however,  the  wave  front  is  typically  not  pure  TEMq0 .  The 
phase  angle  is  therefore  not  constant,  but  is  specified  in  terms  of  the 
local  field  0(r,z)  by 


lD0,0(0'r)l 

sin®(r,z)  »  sin®  — — - — — 

o  I D( r,z )  I 


sin®^ 
_ o 

I  U(  r,  z )  I 


[11] 


where  Zr  «  n w*/X,  wQ  is  the  1/e  E-field  radius  of  the  TEM00  beam  at  the 
waist,  lu00(0,z)l  represents  the  on-axis  amplitude  of  a  TEMq0  beam 
normalized  for  unity  amplitude  at  the  waist, and  ®Q  is  the  stable-phase 
angle  which  would  be  produced  on  axis  in  the  presence  of  a  TEM^  beam  of 
this  normalization.  The  initial  pump  beam  is  a  TEMQ0  wave  with  unity 
normalization.  As  the  modal  structure  evolves  due  to  the  FEL  interaction, 
the  resonant  phase  evolves  self-consistently  according  to  Eq.  [11]. 
Relatively  high  E-field  values  cause  the  electrons  to  slide  toward  the 
bottom  of  the  bucket  (toward  small  phase  angles  with  high  refraction  and 
low  gain)  while  relatively  small  E-fields  cause  large  phase  angles  and 
high  gain.  This  effect  tends  to  produce  an  equilibrium  value  of 
intracavity  flux. 


In  principle,  the  phase  angle  varies  radially  within  the  e-beam. 

Such  variation  has  been  included,  but  is  unimportant  for  the  cases 
examined  since  the  photon  intensity  is  nearly  constant  across  the  e-beam. 
If  the  phase  angle  reaches  90  degrees,  electron  detrapping,  or  escape  from 
the  ponderomotive  potential  well  and  subsequent  loss  of  energy  exchange, 
will  occur.  We  have  allowed  electron  detrapping  to  occur  as  a  result  of 
radial  variation  in  the  photon  intensity  distribution,  but  this  is  also 
unimportant  for  these  cases.  The  code  treats  detrapping  by  assuming 
straight-line  electron  trajectories  parallel  to  the  wiggler  axis.  This 
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neglects  the  sinusoidal  betatron  orbits  [12]  which  result  fron  transverse 
gradients  in  the  magnetic  field.  The  code  is  nevertheless  applicable  to 
cases  which  have  significant  betatron  motion,  provided  that  there  is 
little  radial  variation  of  photon  intensity  and  phase  across  the  e-beam. 

Detrapping  may  also  occur  across  the  entire  e-beam  as  a  result  of 
mode  beating.  For  example,  it  is  shown  later  that  interference  between 
the  1-0  and  0-0  modes  can  lead  to  axial  intensity  fluctuations  of  roughly 
a  factor  of  two.  If  a  design  phase  angle  of  BQ  -  45  degrees  had  been 
chosen.  Eg.  [11]  shows  that  detrapping  would  occur  on  axis  when  the 
competition  of  the  two  modes  is  included.  Since  0(r,z)  is  calculated 
based  on  the  local  field,  any  detrapping  caused  by  mode  beating  is 
self-consistently  taken  into  account.  The  detrapping  problem  can  be 
avoided  by  choosing  a  smaller  design  value  of  the  phase  angle,  but  at  the 
expense  of  electron  deceleration. 

The  oscillator  mode  structure  may  be  characterized  by  decomposing 
the  wave  front  into  the  axisymmetric  Gaussian-Laguerre  modes  of  the 
resonator.  The  normalized,  orthogonal  modes  may  be  written[9] 

Vo->  ->f  1  lpl2fe)  1  112 

•  exp[-i  {*»  -  <2p  *  1)  fn“  [i]  *  i  (jjy  ]] 

where  w(z)  -  w^l  +  (z/ZR)2,  p  -  0,1,2,...  is  an  integer  index  for  each 
mode,  the  Lp  are  the  associated  Laguerre  polynomials  of  order  p,  and 

2  A  A 

ZR  -  pw0/X  is  set  by  the  mirror  curvature.  The  inodes  U0  0,  0, 

A 

U2  0,...  correspond  to  circular  TEMq0,  TEM1q,  TEM20,...  modes, 
respectively.  As  the  modes  propagate  through  one  round  trip  of  the 
cavity,  they  incur  a  phase  shift  relative  to  a  plane  wave  of 
4(  2p+l  )tan”  *( L/2Zg )  due  to  curvature. 


We  can  describe  any  arbitrary  wave  front  as  an  expansion  of  nodes 


&(r,r)  «  E  »p  Up  0<r,*)  .  [13] 

p-o 

The  complex  coefficients  are  conveniently  found  by  evaluating 

*  -  - 

a  -  I  U(r,z)U  _(r,z)  277rdr  ,  [14] 

y  *  q  y  • v 

where  the  superscript  denotes  a  complex  conjugate.  On  each  pass 
through  the  laser  the  wave-front  mode  structure  is  analyzed  at  the  output 
mirror,  immediately  following  the  mirror  truncation.  This  decomposition 
describes  both  the  internally  reflected  and  the  outcoupled  field.  The 
far-field  distribution  may  be  found  by  passing  the  output  beam  through  a 
collimating  lens,  then  focusing  with  a  spherical  mirror.  A 
Huygens -Fresnel  integral[9]  is  used  for  propagation  to  the  focal  spot, 
yielding  the  far-field  pattern. 

Most  of  the  simulations  presented  here  pertain  to  confocal  cavities, 
that  is  the  Rayleigh  range  ZR  is  equal  to  half  the  cavity  length,  Lc. 
Confocal  systems,  as  well  as  concentric  (LC/ZR  «  00 )  and  planar 
(1>C/Z£  -  0),  can  support  unusual  mode  structure  because  the  relative  phase 
slippage  between  cavity  modes  over  one  round  trip  is  an  integral  multiple 
of  2R .  For  empty  cavities,  relative  phases  of  individual  modes  at  any 
given  transverse  plane  will  repeat  on  each  round  trip.  Consequently,  the 
mode  mixing  due  to  truncation  by  any  apertures  present  has  interesting 
properties.  For  example,  the  gain-free  confocal  system  injected  with  pure 
0-0  mode  develops  almost  no  1-0  mode  by  way  of  mirror  truncation,  but  the 
amount  of  2-0  mode  can  be  appreciable.  The  1-0  mode  is  produced  by 
truncation  of  the  0-0  mode  at  each  mirror,  but  the  1-0  fields  generated  at 
the  two  mirrors  drive  1-0  modes  which  are  1B0  degrees  out  of  phase.  The 
2-0  mode  fields  produced  at  each  mirror  are  in  phase.  Zf  1-0  mode  is 
produced  by  the  gain  media,  another  interesting  effect  occurs.  In 
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propagating  from  one  end  of  the  cavity  to  the  other,  the  phase  of  the  0-0 
■ode  changes  by  v/z  (relative  to  a  plane  wave).  While  that  of  the  1-0  mode 
■lips  by  3fr/2.  The  interference  is  therefore  different  at  each  end  of  the 
cavity,  providing  an  asymmetry  in  the  axial  direction.  Such  an  asymmetry 
can  be  supported  in  the  PEL  since  the  gain  mechanism  is  active  in  one 
direction  only. 

These  effects  are  less  apparent  in  cavities  intermediate  to  the 
special  cases  of  confocal,  concentric,  or  planar.  The  importance  of  any 
deleterious  effects  might  therefore  be  considered  minimal,  but  the  reverse 
may  actually  be  true.  The  low  Fresnel  number  and  small  wiggler  bores 
inherent  to  FELs  tend  to  result  in  unusually  high  energy  loading  on  mirror 
surfaces.  For  high  average  power  systems,  alleviation  of  these  problems 
■ay  require  the  use  of  very  long,  near-concentric  cavities. 

XV.  RESULTS 

4.1  Properties  of  a  Confocal  PEL  Resonator 

The  mode  structure  of  the  confocal  FEL  oscillator  described  in  Table 
I  has  been  investigated  in  detail.  We  will  first  point  out  some  of  the 
basic  features  of  the  mode  structure  of  a  confocal  FEL  oscillator, 
followed  by  an  investigation  of  the  effect  of  varying  some  of  the 
parameters . 

The  conditions  listed  in  Table  Z  describe  a  confocal  stable 
resonator  with  apertures  of  radius  1.8  w,  where  w  is  the  1/e  radius  of  the 
E-field  of  a  TEMqq  mode  at  the  ends  of  the  wiggler. 
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» 
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TABLE  I 


Standard  Conditions  for  PEL  Oscillator 
Mode  Structure  investigation 


Cavity  Length,  L^Z*  2  (confocal) 

Aperture  Size,  a^/w  1.8 

E-Beam  Size,  Z^/Z^  10 

Design  Phase  Angle,  eo  20.7* 

Gain  10%/pass 

Output  Coupling  9 . 1% 

Aperture  Truncation  Length,  T/aw  0.055 


This  mirror  size  provides  a  modest  degree  of  mode  selection.  Higher  order 
modes  are  attenuated  since  they  carry  more  energy  at  larger  radii.  The 
electron  beam  size  is  specified  by  Zj^/zJ^  -  io,  where  we  define 
zj^  -  jrrJj/X  to  be  the  Rayleigh  range  based  on  the  l/e  radius  of  the 
Gaussian  e-beam.  We  choose  the  electron  beam  size  to  be  significantly 
smaller  than  the  pump  beam  (*eb/wQ  -  0.316),  so  the  e-beam  is  essentially 
pumped  by  a  uniform  intensity  wave  front.  The  peak  (on-axis)  electron 
density  is  chosen  to  give  a  gain  approximately  equal  to  the  9.1  percent 
output  coupling  (reflectivity  of  0.909).  Small  additional  cavity  loss 
results  from  clipping  at  the  apertures. 

The  distribution  of  Gaussian-Laguerre  modes  at  the  output  coupler 
are  shown  as  a  function  of  round-trip  number  in  Pig.  2. 

The  wave  front  is  initiated  as  a  TEMq0  mode,  but  settles  down  within  30 
passes  to  a  steady  mode  distribution  containing  a  significant  contribution 
of  1-0  mode.  All  higher  modes  are  strongly  suppressed  by  the  finite 
apertue  radii.  The  high  content  of  1-0  mode  is  due  to  the  PEL 
interaction,  rather  than  by  truncation  at  the  apertures.  The  fraction  of 

1- 0  mode  power  in  the  absence  of  gain  is  only  2.6  x  io-*.  The  higher 
order  mode  with  the  highest  energy  content  in  the  bare  resonator  is  the 

2- 0  mode;  it  constitutes  only  0.27  percent  of  the  total  power. 
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Figure  2.  Mode  structure  evolution  for  conditions  of  Table  I. 

The  electron  trapping  fraction  and  bean  power  are  ahown  as  a 
function  of  round-trip  nunber  in  Fig.  3. 


Figure  3.  Beam  power  and  fraction  of  electrons  trapped  as  a  function 
of  round-trip  number  for  conditions  of  Table  I. 


As  the  aode  structure  evolves,  the  round-trip  losses  exceed  the  gain 
slightly  on  each  pass,  allowing  the  power  level  to  decay  about 
25  percent.  Once  a  steady  mode  structure  has  developed,  the  power  level 
cosies  to  steady  state  as  the  phase  angle  adjusts  to  allow  the  gain  to 
match  the  round-trip  losses.  The  electron  beast  is  trapped  essentially 
over  its  full  radius. 

The  steady-state  phase  angle  of  the  PEL  interaction  is  shown  as  a 
function  of  axial  and  radial  position  in  Pig.  4. 


Figure  4 .  Steady-state  phase  angle  for  FEL  interaction  as  a  function  of 
position.  Conditions  given  in  Table  Z. 

While  the  on-axis  phase  angle  is  initially  20.7  degrees  for  the  0-0  node 
pump  beam,  the  steady  node  distribution  results  in  on-axis  phase  angles 
which  vary  between  15  and  32  degrees.  The  largest  on-axis  phase  angle,  or 
lowest  intensity,  occurs  at  the  back  mirror.  The  radial  phase  angle 
profiles  have  a  rather  broad,  flat  distribution,  with  a  sharp  transition 
to  large  angles  and  detrapping  at  large  radius.  The  broad  distribution  is 
consistent  with  the  e-beam  being  pumped  by  an  essentially  uniforn 
intensity  wave  front.  The  narrowest  profile  occurs  at  the  midplane 
between  the  mirrors  (at  the  photon  beam  waist)  where  detrapping  occurs  at 
r  ■  2.4  r_h .  Downstream  from  this  point,  all  electrons  at  this  radius  or 


Figure  5.  Intracavity  intensity  distributions  of  steady  mode  structure 

on  (a)  forward  pass,  and  (b)  return  pass  for  conditions  of  Table  I. 

The  intensity  distribution  I(r,z)  for  the  forward  pass  is  shown  in 
Fig.  5(a)  anJ  the  backward  pass  in  5(b).  The  intensity  is  normalised  in 
such  a  way  that  the  on-axis  intensity  for  the  injected  0-0  mode  is  1  at 
the  midplane  (z/ZR  »  0)  and  0.5  at  the  mirrors  (z/Zg  •  ±1).  The  sharp 
variation  of  intensity  near  the  axis  just  after  a  mirror  reflection  is  due 
to  the  diffraction  pattern  of  the  aperture.  The  steady  distribution 
develops  in  such  a  way  that  the  central  on-axis  intensity  is  higher  than 
that  of  the  individual  0-0  mode  pump  beam  on  forward  passes  and  lower  on 
backward  passes.  The  forward  and  backward  going  on-axis  intensities  are 
compared  with  the  0-0  mode  in  Fig.  6. 


2 


Figure  6.  Comparison  of  on-axis  intensity  distributions  of  steady  mode 
structure  and  injected  0-0  mode. 


The  on-axis  intensity  on  the  output  mirror  is  over  two  times  that  of  the 

back  mirror.  This  difference  cannot  be  accounted  for  by  gain,  which  is  • 

only  10  percent,  but  results  from  interference  between  the  dominant  0-0 

and  1-0  modes.  The  relative  phases  of  the  two  modes  are  such  that  they 

constructively  interfere  on-axis  on  forward  passes  and  at  the  output 

mirror  while  destructively  interfering  on  backward  passes  and  at  the  back  • 

mirror . 

The  radial  structure  of  the  output  beam  is  shown  by  the  intensity 
and  phase  profiles  in  Fig.  7.  9 


-•  • 
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RADIUS,  r/a 

Figure  7.  Intensity  and  phase  profiles  at  output  coupler  for  conditions 
of  Table  I.  Phase  profile  has  spherical  curvature  of  output 
mirror  removed.  Intensity  profiles  are  normalized  for 
equal  power. 


The  conditions  are  those  of  Table  I,  except  that  data  is  included  for 
outcoupling  at  either  end  of  the  cavity.  The  axial  structure  and 
Gaussian-Laguerre  mode  content  of  the  two  cases  is  essentially  identical; 
the  differences  in  the  intensity  and  phase  profiles  may  be  attributed  to 
different  relative  phases  between  the  0-0  and  1-0  modes.  When  the 
spherical  curvature  of  the  output  mirror  is  removed  from  the  phase 
distributions  of  the  two  nearly  Gaussian  output  beams,  each  beam  contains 
some  additional  residual  spherical  curvature.  Evidently  the  beam  at  the 
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•-bean  exit  has  somewhat  nor*  spherical  curvature  than  the  output  mirror 
and  the  bean  at  the  e-bean  entrance  has  somewhat  less  curvature. 

A  useful  measure  of  the  focusability  of  the  beans  is  made  by 
examining  their  amplitude-weighted  rms  phase  deviations  at  the  output 
mirror.  When  all  spherical  curvature  is  removed  from  the  phase  profiles, 
the  residual  phase  aberrations  shown  in  Fig.  8  are  found. 


RADIUS,  r/s 

Figure  8.  Residual  wave-front  distortion  at  output  mirror  with  all 
spherical  curvature  removed. 

All  phase  differences  across  the  aperture  are  less  than  l/io  of  a 
wavelength.  Since  the  aberrations  are  predictable  rather  them  random, 
they  could  be  corrected  by  use  of  properly  designed  optical  components. 
Such  an  effort  is  probably  unnecessary,  however,  since  the  aberrations  are 
so  small.  The  amplitude-weighted  rms  wave-front  distortion[13]  is 
approximately  1/50  of  a  wavelength  at  both  the  e-beam  exit  and  entrance. 
This  phase  deviation  results  in  a  far-field  peak  intensity  loss  of  less 
than  2  percent.  Equivalently,  the  Strehl  ratio  is  greater  than 
98  percent. 


An  Additional  demonstration  of  output  bean  quality  is  the 
calculation  of  power  delivered  to  a  far-field  target.  Conceptually,  such 
a  test  is  nade  using  the  laser  together  with  a  focusing  lens.  An 
appropriate  figure  of  nerit  is  the  power  which  can  be  delivered  within  a 
specified  target  area  using  a  lens  of  specified  aperture  size.  A  nearly 
equivalent  mathematical  figure  of  merit  involves  the  fraction  T)  of  total 
power  Fjot  which  can  be  delivered  within  a  radius  R^  on  the  target 

2  R. 
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Here  f  is  the  distance  from  focusing  lens  to  target,  and  a  is  the  mean 
radius  of  the  beam  as  it  exits  the  focusing  optics. 


a  ■  — - —  J" r  l(r)2rrrdr. 
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Results  are  shorn  in  Fig.  9  as  a  function  of  the  target  radius.  The 
target  radius  has  been  normalized  to  a/f. 


NORMALIZED  TARGET  RADIOS,  R,/(fX/2n*) 

Figure  9.  Fraction  of  power  enclosed  at  focal  plane  as  a  function  cf 
normalized  radius. 


Results  for  a  diffraction-limited  TEMqq  beam  are  included  for  comparison. 
The  difference  in  power  delivered  to  the  target  for  these  two  output  beams 
is  very  small.  The  high  focusability  of  each  beam  is  evident. 

Such  a  conclusion  may  appear  to  contradict  the  mode  decomposition  of 
Fig.  2,  which  shows  about  20  percent  of  the  energy  to  be  in  the  1-0  mode. 
The  apparent  discrepancy  lies  in  the  fact  that  Fig.  2  is  a  decomposition 
based  on  the  basis  set  determined  by  the  mirror  geometry.  The  choice  of 
different  basis  sets,  i.e.,  a  change  of  Rayleigh  range  and  waist  location, 
will  result  in  different  energy  partitions.  We  have  made  an  empirical 
search  of  Rayleigh  range  and  waist  location  parameters  for  the  purpose  of 
finding  the  highest  projection  into  the  0-0  mode.  The  result,  for  the 
case  of  outcoupling  at  the  e-beam  exit,  is  that  96.7  percent  of  the  total 
energy  can  be  projected  into  the  0-0  mode  using  the  best  matched  Rayleigh 
range  and  waist  location.  Physically,  this  means  that  the  intracavity 
beam  is  of  nearly  diffraction-limited  quality,  but  its  curvature  is  not 
matched  to  that  of  the  fundamental  cavity  mode.  The  mismatched  wave  front 
does  not  reflect  back  on  itself,  resulting  in  the  front-to-back 
asymmetry . 

4.2  Parametric  variation 

we  have  seen  that  for  a  nominal  phase  angle,  0O,  of  20.7  degrees, 
approximately  22  percent  of  the  total  energy  of  the  stable  transverse  mode 
lies  in  the  1-0  cavity  mode.  The  phase  angle  is  a  design  parameter, 
depending  primarily  on  the  wiggler  taper,  cavity  Q,  and  e-beam  current. 
Analysis  of  the  mode  structure  variation  as  a  function  of  BQ  is  of 
interest,  since  a  range  of  design  values  will  be  achievable  in  practice. 
This  has  been  investigated  with  two  simulations  in  which  phase  angles  of 
90  degrees  and  0  degrees  were  used.  These  two  limiting  cases  correspond 
to  pure  gain  and  pure  refraction,  respectively.  In  both  cases,  the  phase 
angle  was  fixed  everywhere  in  r  and  z  so  that  there  is  no  possibility  of 
electron  detrapping.  All  other  parameters  are  those  of  Table  I.  The 
results  are  shown  in  Fig.  10. 
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Figure  10.  Mode  structure  evolution  for  various  FEL  interaction  phase 
angles  with  fixed  trapped  electron  density.  Phase  angle 
is  everywhere  uniform  with  no  electron  detrapping. 

ZR/ZRb  "  10 '  av/w  "  1*80'  Lc/ZR  *  2’ 


When  compared  to  the  earlier  results,  the  pure  refraction  case  develops  a 
large  fraction  of  higher  order  modes.  In  fact,  the  round-trip  diffraction 
losses,  i.e.,  power  falling  outside  the  mirror  radius,  approach 
10  percent.  Evidently,  refraction  is  a  much  strc  ger  mechanism  for  beam 
spreading  than  gain.  While  large  phase  angles  appear  to  be  desirable  for 
minimizing  the  higher  mode  content,  there  are  fundamental  limitations  on 
how  large  a  design  phase  angle  can  be  used.  For  example,  as  discussed  in 
Sec.  Ill,  larger  design  phase  angles  are  more  susceptible  to  detrapping 
due  to  mode  beating. 
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This  variation  in  phase  angle  was  Bade  with  a  fixed  trapped  electron 
density.  It  could  he  achieved  in  practice  by  variation  of  the  cavity  Q, 
thereby  changing  the  cavity  flux  level  and  rotating  the  locally  generated 
electron  E-field  phasor  with  respect  to  the  applied  E-field  phasor.  A 
related  variation  is  the  change  in  length  of  the  electron  E-field  phasor 
with  the  angle  fixed. 

Results  are  shown  in  Pig.  11. 
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Zn  this  case,  use  of  low  output  coupling  results  in  a  general  reduction  in 
higher  order  mode  content,  since  both  gain  and  refraction  are  reduced.  As 
the  gain  is  increased  above  10  percent,  the  content  of  1-0  node  appears  to 
reach  a  plateau.  Apparently,  there  is  a  mechanism  for  either  increasing 
the  losses  of  1-0  node  or  decreasing  its  production,  causing  the  1-0  node 
content  to  saturate  at  a  level  independent  of  the  length  of  the  electron 
E-field  phasor. 

The  effect  of  changing  the  electron-beam  size  with  fixed  gain  is 
shown  in  Pig.  12. 


Figure  12.  Mode  structure  for  various  e-beam  sizes  with  fixed  gain. 

10  percent,  9  percent  output  coupling,  6  *  20.7° 
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As  the  e-beam  diameter  is  varied,  the  trapped  electron  density  is  changed 
so  as  to  keep  the  gain  constant.  The  e-beam  size  is  not  increased  to  the 
point  that  detrapping  of  the  wings  of  the  e-beam  density  distribution 
becomes  significant.  Increasing  the  size  of  the  e-beam  decreases  the 
content  of  higher  modes.  This  is  consistent  with  the  notion  that  a  small 
diameter  gain  media  leads  to  high  diffraction  angles  and  hence  high  order 
modes.  However,  in  this  case,  the  e-beam  constitutes  an  antenna  of 
sufficient  length  that  the  angular  spread  of  its  emission  is  nearly  length 


C-23 


dominated  and  independent  of  transverse  size.  Accordingly,  we  expect  that 
reduction  of  the  e-beam  size  to  zero  will  not  produce  mode  structure 
dramatically  different  from  that  of  the  Z^/zJ*5  -  20  curve  of  the  figure. 

The  steady-state  mode  structure  is  also  a  function  of  the  aperture 
size.  As  shown  in  rig.  13,  the  use  of  small  apertures  strongly  suppresses 
higher  order  modes. 


HOUND  1W*  WWW 

Figure  13.  Mode  structure  for  several  aperture  sizes.  60  *  20.7®, 

10  percent  gain,  9  percent  output  coupling,  ZR/Zfb  *  10, 

VS  "  2‘ 

The  diffraction  losses,  defined  as  the  fraction  cf  the  beam  power  falling 
outside  the  mirror  radius,  actually  decrease  for  the  smaller  mirror  size. 
The  round-trip  fractional  diffraction  loss  is  0.0310  for  a^/w  -  1.80  and 
0.0221  for  Sy/w  -  1.55.  Apparently  there  is  a  finite  aperture  size  for 
which  the  diffraction  losses  are  minimized.  For  the  infinite  mirror  size 
there  is  no  mode  selection,  since  the  output  coupling  is  mode 
independent.  Zn  this  case,  it  is  interesting  to  note  that  the  relative 
content  of  0-0  and  1-0  modes  oscillates. 

The  last  variation  considered  is  that  of  cavity  length.  As 
discussed  in  Sec.  ZZZ,  mode  beating  in  confocal  and  concentric  cavities 


has  special  properties  since  the  round-trip  phase  slippage  between  Bodes 
is  an  integer  aultiple  of  2v.  This  is  seen  in  Fig.  14  where  confocal  and 
concentric  cavities  support  large  fractions  of  TEMi0  node,  but  for 
intermediate  cavity  lengths  the  fraction  of  1-0  node  decreases 
dramatically. 
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Figure  14.  Mode  structure  for  various  cavity  lengths,  6o  *  20.7°, 

ZR/Z|b  *  10,  aw/w  ■  1.80,  a^/w  ■  3,  gain  ■  10  percent, 

9  percent  output  coupling. 

It  was  shown  previously  (Fig.  6)  that  the  presence  of  higher  order  modes 
has  unusual  effects  on  the  intensity  distribution  within  the  wiggler.  In 
the  case  of  the  concentric  cavity,  the  optical  bean,  when  compared  to  the 
fundamental  cavity  mode,  is  relatively  broad  and  low  intensity  where  it 
enters  the  wiggler.  Before  exiting  the  wiggler,  the  beam  then  focuses  to 
a  relatively  narrow,  intense  waist.  This  is  analogous  to  the  behavior 
seen  on  forward  passes  in  the  confocal  cavity  (see  Fig.  6).  But  in 
contrast  to  the  confocal  results,  the  spot  sises  on  the  two  mirrors  and 
the  Intensity  distributions  on  forward  and  reverse  passes  are  very  similar 
in  the  concentric  cavity. 

The  lack  of  higher  order  mode  content  in  cavities  of  intermediate 
length  results  from  the  phase  of  the  1-0  mode  being  determined  by  the 


ponder omotive  potential,  which  is  formed  by  the  dominant  0-0  mode.  Thus, 
the  1-0  phase  is  tied  to  the  0-0  mode,  and  the  1-0  wave  produced  on  a 
given  round  trip  is  not  in  phase  with  1-0  waves  produced  on  earlier  or 
subsequent  round  trips  in  a  cavity  of  intermediate  length.  A  further 
manifestation  of  the  varying  phase  relationship  between  modes  is  that  the 
intensity  distribution  on  the  mirrors,  and  hence,  the  cavity  losses,  are 
variable  from  one  round  trip  to  another.  This  results  in  the  ripple  seen 
in  Pig.  14. 

Nevertheless,  a  steady  mode  structure  develops  which  has  a  constant 
relative  phase  between  any  two  modes  at  any  given  location  on  any  pass. 
This  is  because  the  production  of  higher  order  modes,  either  by  mirror 
truncation  or  PEL  gain  and  refraction,  always  occurs  at  a  constant  phase 
angle  relative  to  the  existing  pump  beam  (i.e.,  coherently).  The  higher 
order  mode  which  is  produced  (say  1-0  mode,  for  example)  will  continue  to 
undergo  phase  slippage  relative  to  the  0-0  mode  as  it  propagates  around 
the  cavity.  However,  it  is  also  gradually  attenuated  and  during  its 
lifetime  it  has  a  fixed  phase  relative  to  the  0-0  mode  at  any  given 
location . 

Of  the  three  cavity  lengths  considered,  the  confocal  and  concentric 
cavities  have  the  highest  content  of  higher  order  modes.  The  round-trip 
diffraction  losses  are  also  highest  for  these  cavities,  as  indicated  in 
Table  II. 


Table  ZZ 

Round-Trip  Aperture  Losses  for 
various  Cavity  Lengths 

Cavity  Length,  L  /Z 

_ C  R 

-  200 

4.8 
2 


Round-Trip  Fractional 
Diffraction  Loss 

0.036 

0.009 

0.031 
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these  losses  pertain  to  energy  passing  beyond  the  edges  of  the  mirrors  and 
energy  lost  due  to  aperturing  by  the  viggler  magnet  or  other  elements  in 
the  beam  line.  For  the  confocal  and  concentric  cases,  the  loss  is  about 
one-third  of  the  output  coupling,  representing  relatively  inefficient 
energy  extraction. 

V.  SUMMIT  MID  CONCLUSIONS 

The  transverse  structure  of  the  tapered-wiggler  PEL  optical  field 
has  been  analysed  by  numerical  solution  of  a  paraxial  wave  equation.  The 
optical  cavity  is  initially  injected  from  an  external  source  at  its  full 
saturated  intensity  and  the  subsequent  transverse  structure  evolution  is 
of  interest.  Unusual  features  of  the  PEL  geometry,  as  compared  to 
conventional  lasers,  include  the  narrowness,  or  more  precisely  the  low 
Presnel  number  of  the  gain  media,  the  small  diameter  of  the  gain  media 
with  respect  to  that  of  the  optical  field,  and  the  lack  of  gain  on  return 
passes.  These  factors  tend  to  produce  a  steady-state  mode  structure 
different  from  the  pure  TEMq.  injected  wave. 

The  higher  order  mode  content  is  especially  evident  in  confocal  and 
near-concentric  optical  cavities,  i.e.,  those  with  mirrors  separated  by 
one  or  two  times  their  radius  of  curvature.  Por  such  cavities,  the 
round-trip  phase  shift  between  nodes  is  an  integer  multiple  of  zv ,  a 
special  situation  resulting  in  constructive  addition  of  higher  order  modes 
produced  on  different  round-trips.  The  presence  of  higher  order  modes  can 
significantly  affect  the  photon  beam  shape,  providing  unusual  effects. 

Por  example,  in  the  confocal  cavity  case,  the  beam  exhibits  different 
radial  structure  on  forward  and  reverse  passes,  as  well  as  different  spot 
sixes  on  the  front  and  rear  mirrors.  Practical  PEL'S  operating  at 
substantial  average  power  levels  are  likely  to  require  near-concentric 
cavities  in  order  to  provide  sufficiently  large  spot  sixes  at  the 
mirrors.  In  this  case,  2»r  shift  occurs  between  the  dominant  TEMq0  and 


TEM10  aodti  in  on*  pass,  rather  than  a  round  trip,  so  that  spot  sizes  on 
the  front  and  rear  airrors  and  forward  and  reverse  propagating  waves 
differ  fro*  those  of  the  fundamental  node  but  are  comparable  to  each 
other. 

The  quality  of  the  output  beam  is  excellent  for  the  cases  studied, 
being  of  nearly  diffraction-limited  quality.  This  result  is  somewhat 
surprising  in  view  of  the  TEM10  mode  content  which  complicates  the 
intracavity  structure.  This  apparent  contradiction  may  be  understood  by 
noting  that  the  intracavity  wave-front  curvature  is  not  matched  to  that  of 
the  fundamental  mode  of  the  cavity,  while  the  intracavity  beam  is 
essentially  diffraction  limited. 

Additional  findings  involve  the  choice  of  synchronous  phase  angle,  a 
parameter  chosen,  within  limits,  at  the  discretion  of  system  designers. 

It  is  found  that  for  a  given  electron  density,  a  synchronous  phase  angle 
chosen  to  maximize  gain  results  in  much  less  TEM10  (or  higher  order)  mode 
production  than  does  a  phase  angle  chosen  to  maximize  the  phase  shift.  In 
the  former  case,  the  e-beam  acts  as  an  extended  antenna  which  produces  a 
field  much  like  that  already  in  the  cavity,  while  in  the  latter  case,  the 
e-beam  acts  like  a  series  of  focusing  lenses. 

A  useful  application  of  the  techniques  described  would  be  analysis 
of  the  effects  of  hole  couling  as  a  means  of  output  coupling  in  cavities 
compatible  with  high  average  power.  Extension  of  these  techniques  to 
include  the  unbunched  nature  of  electrons  entering  the  viggler  and 
subsequent  oscillation  of  the  phase  angle  in  the  ponderomotive  potential 
well,  would  be  useful.  A  related  effect  of  interest  is  the  slippage  of 
photons  relative  to  electrons  when  short  pulses  are  considered.  In  this 
case,  the  transverse  structure  is  modified  due  to  reduction  of  the 
effective  interaction  length. 
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Abstract 

Simple  scaling  laws  for  free-electron  laser  oscillator  mirror  alignment  tolerances  are 
given  based  on  geometric  optics.  The  effect  of  geometric  walk-off  and  diffractive  scraping 
are  discussed.  In  cavities  that  are  nearly  concentric,  the  effect  of  diffraction  can  relax 
the  alignment  tolerances. 


Introduction 


Free-e lec-iron  laser  (FEL)  oscillator  cavities  are  typically  going  to  be  near-concentric; 
the  centers  of  curvature  of  the  cavity  endmirrors  will  nearly  overlap.  As  will  be  explained, 
near-concentric  cavities  are  alignment-sensitive.  In  designing  a  FEL  cavity,  an  optical 
waist  extending  several  meters  or  more  through  the  wiggler  may  be  desirable  for  good  gain 
and  energy  extraction,  giving  the  beam  a  low  far-field  divergence.  The  cavity  must  be  long 
enough  to  avoid  mirror  damage.  For  a  simple  two-mirror  symmetric  cavity,  the  choice  of 
optical  waist  parameters  and  cavity  length  fixes  the  required  radius  of  curvature  of  the 
end  mirrors.  The  alignment  sensitivity  of  this  cavity  is  found  by  determining  when  mis¬ 
alignments  degrade  performance.  Degradation  may  be  due  to  loss  of  electron-beam/optical- 
beam  overlap,  off-axis  optical  aberration,  or  losses  at  apertures  defined  by  the  wiggler 
bore,  the  mirror  size,  or  a  scraping  outcoupler.  In  the  analysis  to  follow,  the  scaling  of 
alignment  sensitivity  with  FFL  parameters  will  be  found  based  on  the  effects  of  electron- 
bear  /optical-beam  overlap  or  losses  at  apertures. 

Length  requirements  and  cavity  design 

The  cavity  design  for  a  FEL  oscillator  is  constrained  by  the  choices  of  optical  power, 
wavelength,  and  the  efficiency  of  the  laser.  Large  electron-beam  energy  extraction  is 
needed  for  high  efficiency.  For  large  extraction,  high  peak  optical  intensities  are  needed 
at  a  long  optical  waist  in  the  wiggler  for  the  duration  of  the  electron-beam  micropulse. 

The  length  of  the  waist  region  is  described  in  terms  of  the  Rayleigh  range,  Zp,  the  distance 
over  which  a  Gaussian  beam  radius  grows  by  a  factor  of  *T.  A  Rayleigh  range  chosen  to 
optimize  laser  performance  is  typically  from  1/2  to  1/4  the  wiggler  length.  Larger  wigglers 
tend  to  give  superior  laser  performance  but  at  increased  wiggler  cost,  at  possibly  increased 
electron-beam  steering  difficulty  in  the  wiggler,  and  as  shown  later,  at  increased  cavity 
length  to  avoid  mirror  damage.  Rayleigh  ranges,  in  present  experiments,  vary  from  0.1  to 
lm,  and  may  be  considerably  larger  in  advanced  oscillator  designs.  Choosing  the  Rayleigh 
range  and  wavelength  determines  the  shape  of  an  unaberrated  freely-propagating  Gaussian  beam 
diverging  from  a  waist. 

The  minimum  cavity  length,  L,  is  computed  below  based  on  the  requirement  that  the  optical 
beam  radius  at  the  endmirror  be  large  enough  to  keep  the  optical  flux  below  the  damage  or 
degradation  limits  of  the  mirror.  A  simple  two-mirror  cavity  will  be  considered,  as  shown 
in  Figure  1.  Mirror  alignment  tolerances  will  be  derived  subsequently  using  this  analysis, 
and  the  tolerances  will  be  shown  to  hold  even  if  the  cavity  is  shortened  using  a  high  flux 
intracavity  beam  expander,  such  as  with  glancing  incidence  optics.  To  avoid  damage  or 
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Figure  1 


FEL  cavities  are  near-concentric  and  alignment  sensitive. 


degradation  of  the  endmirrors,  which  are  assumed  to  be  used  at  normal  incidence,  the  beam 
radius,  u,  at  the  endmirror  must  often  be  far  larger  than  at  the  waist.  The  constraint  is 
that  the  flux  or  power  density  be  less  than  the  damage  or  degradation  limit 


wu a  >  2I7®d, 
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where  T  is  either  the  average  power  or  single-pulse  integrated  energy  exposure  within  the 
cavity,  whichever  is  damage  limiting;  is  the  damage  limit  for  either  power  density 
(W/cm* )  or  for  single-pulse  integrated  energy  density  (J/cm2),  whichever  is  damage  limiting. 
At  a  distance  L/2  (half  the  cavity  length)  from  the  centered  beam  waist,  the  beam  radius 

is(2) 
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for  L  >>  Zr.  The  wavelength,  X,  Rayleigh  range,  Zr,  and  radius  at  the  waist,  u0, 
are  related  by<2) 
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Using  Equations  1,  2,  and  3,  the  minimal  cavity  length  set  by  damage 
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The  result  of  the  previous  analysis  is  that  a  cavity  may  tend  to  be  very  long.  For 
example,  consider  a  pulsed  oscillator  with  30  ps  micropulses,  and  10*  micropulses  in  a 
macropulse.  If  the  circulating  peak  power  is  S  GW,  Zr  is  lm,  X  is  0.5pm,  and  the  macropulse 
damage  limit  on  a  coating  is  taken  as  10  J/cm*  (which  may  be  conservative  for  very  long 
pulses  but  is  above  typical  for  short  pulse  limits)  then  the  total  flux,  T ,  is  10*  x  30  ps  x 
5  GW  =  1.5  kJ  and  the  required  cavity  length  would  be  about  500m.  This  cavity  length  is 
many  Rayleigh  ranges  long,  and  average  power  requirements  could  lengthen  it  further. 


The  radius  of  curvature  of  the  mirrors  is  designed  to  match  the  optical  wavefront  curva¬ 
ture  for  a  self-replicating  Gaussian  wavefront. (2)  For  a  free-electron  laser,  this  is  a 
good  approximation  but  is  not  exactly  true. (3)  The  radius  of  curvature,  R,  of  the  wavefront 
at  the  cavity  endmirror  is  given  by 

R  "  (§1  (X  *  HTTP']’  (5> 

When  Zr  is  small  compared  to  the  cavity  length 


ZR  «  L/2,  <6> 


then  Equation  5  reduces  to 


R  ■  L/2. 


(7) 


Alignment  tolerance  scaling  using  geometric  optics 

Near-concentric  cavities  are  highly  alignment  sensitive  because  their  geometric  axis, 
defined  by  the  line  joining  the  centers  of  curvature  of  the  two  endmirrors,  undergoes  highly 
leveraged  tilting  when  either  endmirror  is  tilted.  This  is  shown  in  Figure  2  where  the 
cavity  mirror  to  the  right  has  been  tilted  by  an  angle  6,  and  the  optical  axis  is  tilted  by 
an  angle  ♦  .  The  optical  axis  shown  is  almost  cut  by  apertures  at  the  ends  of  the  wiggler 
and  is  misaligned.  If  the  Rayleigh  range  is  much  less  than  the  cavity  length,  the  centers 
of  curvature  of  the  mirrors  are  spaced  by  far  less  than  the  wiggler  length,  and  the  optical 
axis  can  be  considered  to  pivot  about  the  middle  of  the  wiggler.  The  alignment  tolerances 
on  the  endmirrors  can  be  defined  by  requiring  the  optical  axis  stay  within  some  fraction  of 
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Figure  2.  Near-concentric  resonators  are  alignment  sensitive. 


“o  <u0'3>  *s  use^  here)  of  the  Higgler  axis  at  the  Higgler  ends,  thus  preserving  electron 
beam  overlap  and  eliminating  losses  at  the  apertures.  The  tolerance  is 
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nhere  Jw  is  the  length  of  the  Higgler  and  tu0  is  again  the  beam  radius  at  the  naist.  From 
Equations  S  and  7,  it  is  apparent  that  the  mirror  centers  of  curvature  are  displaced  from 
the  cavity  center  by  2Zp/R,  so  that  the  mirror  tilt  (8)  and  the  optical  axis  tilt  (4)  are 
related  by 
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Assuming  the  Higgler  length  is  about  tno  Rayleigh  ranges: 

iw  -  •  2Z K,  (10 

Then  a  tolerance  on  8  can  be  found  in  terms  of  Zp,  1,  ®D,  an(*  r  for  a  cavity  nhere  the 
length  (L)  is  the  minimum  allowed  to  avoid  mirror  damage.  Using  Equations  4  and  7  to  elim¬ 
inate  P  and  L  from  Equation  9,  and  Equations  3  and  10  to  eliminate  u>0  and  tw  from 
Equation  8,  the  tolerance  on  mirror  tilt  is: 
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This  is  the  main  scaling  equation  for  alignment  tolerances  based  on  geometric  optics.  The 
mirror  alignment  sensitivity  is  seen  to  be  tightened  at  short  navelengths,  short  Rayleigh 
ranges,  and  at  higher  poners  or  energies. 

The  tolerance  on  alignment  is  derived  using  geometric,  or  ray,  optics.  The  tolerances 
do  not  change  if  a  telescope  is  inserted  betneen  the  Higgler  and  the  endmirror  in  order  to 
shorten  the  cavity.  The  same  endmirrors  nould  be  used,  having  the  same  radius  of  curvature 
and  area.  The  alignment  sensitivity  depends  only  on  the  alloned  tilt  for  the  optical  axis, 
the  relative  location  of  the  centers  of  curvature  of  the  tno  mirrors,  and  the  degree  to 
which  the  mirror  centers  of  curvature  are  displaced  nhen  the  mirrors  are  tilted.  To  first 
order,  the  alloned  tilt  of  the  optical  axis  is  unchanged  by  use  of  a  telescope,  and  the 
distance  between  the  centers  of  curvature  of  the  two  mirrors  is  also  unchanged.  As  shown 
in  Figure  3,  the  lateral  displacement  of  the  center  of  curvature  of  a  mirror,  when  that 
mirror  is  tilted,  is  also  unchanged  by  a  beam  expander  because  a  magnification,  m,  in  image 
size  corresponds  to  demagnification,  1/m,  in  angle.  The  lateral  displacement  is  given  by 
the  radius  of  curvature  multiplied  by  the  angle  through  which  it  is  tilted.  As  the  cavity 
length  is  demagnified,  the  angle  is  magnified  and  the  displacement  does  not  change  with 
magnification. 
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Figure  3.  Intracavity  telescopes  do  not  affect  the  first  order 
alignment  tolerances  because  a  given  tilt,  6 ,  on  a 
mirror  with  radius  of  curvature  R  results  in  the  same 
displacement  of  the  mirror  center  of  curvature  (6P) 
with  or  without  a  telescope. 


Diffractive  effects 

For  the  example  cited  earlier  of  Zp  ■  100  cm,  X  ■  0.5  x  lO”11  cm,  $D  ■  lo  J/cir*,  and 
T  »  1.5  kj,  the  alignment  tolerance  using  Equation  11  is  about  4  nrad.  The  diffraction 
limit  of  the  mirror  is  roughly  equal  to  the  wavelength  of  the  light  divided  by  the  mirror 
diameter.  The  diameter  is  about  four  times  the  beam  radius  at  the  mirror  and  using 
Equation  1,  the  radius  is  about  11  cm.  The  diffraction  limit  on  pointing  is  about  1  urad, 
over  two  orders  of  magnitude  larger  than  the  alignment  tolerance.  To  sense  such  a  tight 
tolerance,  the  beam  must  circulate  many  times  in  the  cavity. 

If  the  optical  beam  actually  falls  on  the  geometric  optical-axis,  then  the  alignment 
tolerance  given  by  Equation  11  is  valid.  However,  if  a  ray  aligned  to  the  wiggler  axis,  but 
misaligned  to  the  cavity  axis  takes  many  cavity  round  trips  before  seeing  the  losses  assoc¬ 
iated  with  the  misalignment,  then  diffractive  effects  may  cause  the  true  optical  axis, 
defined  by  the  intensity  centroid,  to  be  displaced  from  the  geometric  axis.  To  calculate 
how  long  it  takes  an  off-axis  ray  to  locate  the  geometric  axis,  consider  the  situation 
shown  in  Figure  4.  An  off-axis  ray  is  circulating  paraxially  in  a  near-concentric  bare 
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Figure  4.  An  off-axis  ray  propagating  paraxially  in  a 
near-concentric  cavity. 


cavity.  The  mirror  radius  of  curvature,  R,  is  related  to  its  focal  length,  f,  simply  by 

R  -  2f.  ( 


And  for  a  near-concentric  cavity 


where 


*  “  XT  <<  1- 


Again  L  is  the  cavity  length.  An  off-axis  ray  in  the  cavity  can  be  described  by  a 
vector,  r, 


D-4 


r 


(15) 


where  x  is  the  transverse  position  of  the  ray  and  *'  is  the  transverse  angle.  Propagating 
this  ray  paraxially  can  be  modeled  by  operating  on  the  ray  with  the  matrix,  M,  for  a  cavity 
round  trip. 
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By  examining  the  action  of  this  matrix  on  eigenvectors,  it  can  be  found  that  in  the  limit 
of  large  n 
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where  I  is  the  indentity  matrix.  The  number  of  cavity  round  trips,  n,  needed  for  a  ray  to 
self-replicate  is  therefore  given  by 

ift  n  •  2* i  (It 


n  •  1.5L/Z. 


The  ray  can  be  said  to  find  the  axis  in  roughly  half  this  number  of  passes.  Geometrically 
it  takes  L/ZR  round  trips  for  a  ray  to  find  the  axis.  Diffractive  beam  steering  may  be 
important  if  L/ZR  is  very  large  because  many  diffractive  cuts  occur  prior  to  large 
intensity  loss.  These  effects  were  studied  using  a  diffractive  analysis  of  wavefront 
propagation  in  a  near-concentric  cavity. 

A  model  of  wavefront  propagation  in  a  near-concentric  gain-free  FEL  cavity  was  used  to 
determine  alignment  sensitivity.  In  the  test  case  shown,  the  cavity  is  500  Rayleigh  ranges 
long,  the  wiggler  3  Rayleigh  ranges  long,  and  the  wiggler  bore  4u  at  the  exit.  The 
Rayleigh  range  is  lm.  Optical  propagation  is  modeled  using  a  fast  Fournier  transform 
expansion  of  the  wavefront  into  a  series  of  plane  waves  traveling  at  different  spatial 
angles.  A  TEM00  wave  is  injected  into  the  cavity  and  one  mirror  is  suddenly  tilted. 

The  diffractive  loss  due  to  clipping  at  each  end  of  the  wiggler  is  calculated  and  evolution 
to  a  steady  mode  structure  is  observed.  For  the  best  test  case  shown  in  Figures  5  and  6, 
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Figure  5. 
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Geometric  picture  of  the  case 
studied  to  determine  diffrac¬ 
tive  effects.  The  cavity  was 
misaligned  so  that  the  optical 
axis  was  upon  the  aperture  and 
losses  were  near  100  percent. 
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Fiaure  6.  Optical  lasers  per  round 

trip  for  the  case  shown  in 
Figure  3. 


AS  was  55  nr,  enouqh  to  put  the  geometrical  optical  axis  onto  the  aperture  at  the  wiggler, 
and  create  “100  percent  loss  per  pass.  The  losses  computed  are  far  lower,  equilibrating  at 
about  5  percent  per  round  trip.  Consequently,  the  alignment  tolerance  defined  by  the  geo¬ 
metric  analysis  is  far  too  tight.  This  is  true  in  just  the  cases  when  L/Zr  is  very  large, 
perhaps  10*  or  greater.  Even  so,  the  tolerances  are  still  far  tighter  than  the  diffraction 
limited  pointing  tolerance.  This  is  possible  because  the  laser  resonator,  like  an  etalon, 
allows  the  optical  beam  to  sample  the  mirror  surface  many  times  over  many  round  trips,  and 
more  finely  sense  misalignment.  Diffractive  effects  are  caused  by  the  beam  being  cut  by 
apertures  at  the  edge  of  the  beam  where  the  fields  are  lower,  and  the  intensity  which  is 
proportional  to  the  square  of  the  field,  is  lower  still.  The  cutting  can  steer  the  beam 
without  severe  losses.  The  optical  axis  is  then  not  located  where  the  geometric  analysis 
would  predict. 


Conclusions 


The  alignment  tolerances  for  a  FEL  scillator  built  to  the  minimal  length  allowed  by 
optical  damage  scale  roughly  as  Z#1  X (*d/D  where  Zr  is  the  Rayleigh  range,  X  is  the 
wavelength  and  *p/r  is  proportional  to  the  reciprocal  of  the  mirror  area  required  to  avoid 
optical  damage  due  to  energy  or  average  power  effects.  Short  wavelength,  high  power  FFL 
oscillators  will  have  the  most  demanding  alignment  requirements.  In  near-concentric  FEL 
cavities  which  are  many  Rayleigh  ranges  long,  the  geometrically  derived  tolerances  may  be 
overly  tight  and  diffractive  scraping  at  apertures  will  act  to  keep  the  bear  along  a  lower 
loss  axis.  The  alignment  sensitivity  of  these  cavities  is  always  comparable  to,  or  tighter 
than,  the  diffraction-limited  pointing  accuracy. 
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